MARKOV PROCESSES ON THE PATH SPACE OF THE 
GELFAND-TSETLIN GRAPH AND ON ITS BOUNDARY 



ALEXEI BORODIN AND GRIGORI OLSHANSKI 

Abstract. We construct a four-parameter family of Markov processes on infinite 
Gelfand-Tsetlin schemes that preserve the class of central (Gibbs) measures. Any 
process in the family induces a Feller Markov process on the infinite-dimensional 
boundary of the Gelfand-Tsetlin graph or, equivalently, the space of extreme char- 
acters of the infinite-dimensional unitary group U(oo). The process has a unique 
invariant distribution which arises as the decomposing measure in a natural prob- 
lem of harmonic analysis on U (oo) posed in |Ols03j . As was shown in [BO05a j , this 
measure can also be described as a determinantal point process with a correlation 
kernel expressed through the Gauss hypergeometric function. 
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1. Introduction 

This work is a result of interaction of two circles of ideas. The first one deals with 
a certain class of random growth models in two space dimensions |War07j . |Norl0j . 
BF08+I , |BG09] . IBGR09+I , |BK10] |BorlO+| , while the second one addresses con- 



structing and analyzing stochastic dynamics on spaces of point configurations with 
distinguished invariant measures that are often given by, or closely related to, de- 
terminantal point processes |BU06aj . |BU06bj . |BO09j . jOlslOj . |Olsl0+ . 



Our main result is a construction of a Feller Markov process that preserves the 
so-called zw-measure on the (infinite-dimensional) space Q of extreme characters 
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of the infinite-dimensional unitary group U(oo). The four-parameter family of zw- 
measures arises naturally in a problem of harmonic analysis on U(oo) as the de- 
composing measures for a distinguished family of characters |Ols03] . A zw-measure 
gives rise to a determinantal point process on the real line with two punctures, and 
the corresponding correlation kernel is given in terms of the Gauss hypergeometric 
function |BO05aj . Such point processes degenerate, via suitable limits and/or spe- 
cializations, to essentially all known one- dimensional determinantal processes with 
correlation kernels expressible through classical special functions. 

The problem of constructing a Markov process that preserves a given determinan- 
tal point process with infinite point configurations has been addressed in |Spo87 



|KT10j . | Osa09+| for the sine process, in |KT09] for the Airy process, and in | Olsl0+ | 
for the Whittaker process describing the ^-measures from the harmonic analysis on 
the infinite symmetric group. 

Our approach to constructing the infinite-dimensional stochastic dynamics differs 
from the ones used in previous papers. We employ the fact (of representation the- 
oretic origin) that the probability measures on Q are in one-to-one correspondence 
with central or Gibbs measures on infinite Gelfand-Tsetlin schemes that can also be 
viewed as stepped surfaces or lozenge tilings of a half-plane. The projections of a 
zw;-measure to suitably defined slices of the infinite schemes yield orthogonal polyno- 
mial ensembles with weight functions corresponding to hypergeometric Askey-Lesky 
orthogonal polynomials . 

These orthogonal polynomials are eigenfunctions for a birth and death process 
on Z with quadratic jump rates; a standard argument then shows that the N- 
dimensional Askey-Lesky orthogonal polynomial ensemble is preserved by a Doob's 
^-transform of N independent birth and death processes. 

We further show that the Markov processes on the slices are consistent with re- 
spect to stochastic projections of the iVth slice to the (N— l)st one (these projections 
are uniquely determined by the Gibbs property). This consistency is in no way ob- 
vious, and we do not have a conceptual explanation for it. However, it turns out to 
be essentially sufficient for defining the corresponding Markov process on Q. 

We do a bit more — using a continuous time analog of the general formalism of 
|BF 08+| (which was based on an idea from [DF90] ). we construct a Markov process 
on Gelfand-Tsetlin schemes that preserves the class of central (=Gibbs) measures 
and that induces the same Markov process on Q. 

We now proceed to a more detailed description of our work. 



1.1. Gelfand-Tsetlin graph and its boundary. Following |Wey39| , for iV > 1 
define a signature of length N as an iV-tuple of nonincreasing integers A = (Ai > 
• • • > A at), and denote by GTat the set of all such signatures. Elements of GTn 
parameterize irreducible representations of U(N) or GL(N,C), and they are often 
called highest weights. 
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For A G GUV and v G GT/v+i, we say that A -< v if Uj +i < Xj < Uj for all 
meaningful values of indices. These inequalities are well-known to be equivalent 
to the condition that the restriction of the //-representation of U (N + 1) to U(N) 
contains a A-component. 

Set GT = Utv>i G1V, and equip GT with edges by joining A and v iff A -< v or 
v -< A. This turns GT into a graph that we call the Gelfand-Tsetlin graph. A path 
of length M G {1, 2, . . . } U {oo} in GT is a length M sequence 

A (1) -< A (2) X u) G GT,. 

Equivalently, such a path can be viewed as an array of numbers {Ap } satisfying 

the inequalities aH_^ < A^ < Ap +1 ^; it is also called a Gelfand-Tsetlin scheme. An 
interpretation of paths in GT in terms of lozenge tilings or stepped surfaces can be 
found in the introduction to | BF08+ |. 

The Gelfand-Tsetlin schemes of length N parameterize basis vectors in the Gelfand- 
Tsetlin basis of the irreducible representation of U(N) corresponding to A^, cf. 
[Zhe70]. Denote by Dimjy A the number of such schemes with A^ = A; this is also 
the dimension of the irreducible representation of U(N) corresponding to A. It is 
essentially equal to the Vandermonde determinant in shifted coordinates of A: 

Dim^A) = const^ ] [ (Aj — i — Xj + j). 

l<i<j<N 

A probability measure on infinite paths in GT is called central (or Gibbs) if any 
two finite paths with the same top end are equiprobable, cf. [Ker03| . Let Pn be 
the projection of such a measure to A^ G GTjv. Centrality is easily seen to be 
equivalent to the relation /ijv = /ijv+]A^ +1 , N > 1, where ji^ and /ijv+i are viewed 
as row-vectors with coordinates {hn{.X)}\ & <qt n and {/iAr +1 (z/)} ! , g GT ]V+1 , and 

A ^ +1 = DinMA) i a ^ GTjy, v G GT N+1 , (1.1) 

DimTv(^) 

is the stochastic matrix of cotransition probabilities. There is a one-to-one cor- 
respondence between central measures on GT and characters of U(oo) (equiva- 
lently, equivalence classes of unitary spherical representations of the Gelfand pair 

(U(oo) x [/(oo) diag[/ (oo))), see Pafig j. 

As shown in [Ols03j . see also |Voi76j . |VK82j . jOQ98j . the space of all central 
probability measures is isomorphic to the space of all probability measures on the 
set Q C R^° +2 consisting of the sextuples u = (a+, (3 + , a~, (3~ , 6 + , 5~) G M^° +2 
satisfying the conditions 



a = [a 



>a± >•••>()), /3 ± = (/3f>/3 2 ± >--->0), 5±>0, 

oo 



8=1 
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The set f2 is called the boundary of GT; its points parameterize the extreme 
characters of [/(oo). The map from central measures on GT to measures on Q 
amounts to certain asymptotic relations described in Subsection 19.11 below. 

1.2. ^-measures. Let z, z', w, w' be four complex parameters such that 

(z + k){z' + k) > and (w + k) (V + k) > for any keZ (1.2) 

and 

z + z' + w + w'>-l (1.3) 

(note that (11.21) implies that z + z' and w + w' are real). For N > 1, define a 
probability measure on G1V by (below U = Xi + N — i) 

N 

M z<z<WjW/ \ N (\) = const N |~| {k - lj) 2 Y\_ W z,z>,w, w i(h), (1-4) 

l<i<j<N i=l 

where 



2,2 v ' r(z + - s)r(z / + - x)T(w + 1 + x)T(w' + l + x)' 

We call it the Nth level zw -measure. It is the iV-point orthogonal polynomial 
ensemble with weight W(-), see e.g. |Kon05] and references therein for general 
information on such ensembles. 

One can show that the finite level zw;-measures are consistent: For any N > 1, 
M z>zliW>wl \ N = M z ^^ )Wl \ N+1 A^ +1 . Therefore, the collection {M ZtZ ,^ iW ,\ N } N ^ defines 
a central measure on the paths in GT and a character of U {oo). For z' = z, w' = w, 
this character corresponds to a remarkable substitute for the nonexisting regular 
representation of U(oo), see [Qls03j for details. 

The corresponding measure M ZjZ > w>w i on Q is called the spectral zw-me&sme. If 
co = (a 1 * 1 ,/^,^) G Q is distributed according to M ZiZ ' jWtW / then the random point 
process generated by the coordinates 

\ 2 % 5 2 ' 2 ^« ' 2 » J j=l 

is determinantal, see |BU05aj . [BO05bJ for details. 

1.3. Doob's transforms of A^-fold products of birth and death processes. 

It is not hard to show that the first level zw-measure M Z}Z i )W)V} /\i on GTi = Z is the 
symmetrizing measure for the bilateral birth and death process that from a point 
x G Z jumps to the right with intensity (x — u)(x — u') and jumps to the left with 
intensity (x + v)(x + v'), where (u, u', v, v') = (z, z', w, w'). Denote by V = T> u u i v v i 
the corresponding matrix of transition rates. 
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More generally, we show that the iVth level zw-measure ( 11.4}) is the symmetrizing 
measure for a continuous time Markov chain on GUV with transition rates 

+ £>(V, rijv)l{i i=nij j^ JV }') - d N ■ l x =u (1-5) 

where Zj = Xj + N — j,nj = Uj + N — j, I < j < N, gV is a suitable constant, and 
we take (w, it', t> , v ') = (z + N — 1, z' + N — 1, w, w') in the definition of V. 

Observe that P (7V) can be viewed as a version of Doob's /i-transform of iV copies 
of the Markov chain defined by T> with h(-) = Dinijv( • )■ Note that in our case, 
Dimjy( ■ ) is an eigenfunction of the corresponding matrix of transition rates with a 
nonzero eigenvalue. 

For any N > 1, let (-P/v(t))i>o be the Markov semigroup corresponding to the 
matrix T>^ of transition rates on GTjv (we show that (-P/v(£))t>o is uniquely defined 
and it possesses the Feller property). The key fact that we prove is the consistency 
(or commutativity) relation 

P N+1 {t)A» +1 = A» +1 P N (t), t>0, N>1. 

Although this relation looks natural, we have no a priori reason to expect it to 
hold, and we verify it by a brute force computational argument. 

1.4. Main result. We prove that for any (z, z',w,w') E C 4 subject to (ll.2p - (ll.3l) . 
there exists a unique Markov semigroup (P(t)) t >o on Q that preserves the spectral 
zw-mea.su.ie M z>z > w w i, and whose trace on GTjy coincides with Doob's transforms 
(-P/v(£))t>o introduced above. Moreover, the semigroup (P(i)) t > is Feller (it pre- 
serves Cq(Q), the Banach space of continuous functions vanishing at infinity; note 
that the space Q is locally compact). 

By general theory, see e.g. [EK861 IV.2.7], this means that for any probability 
measure fi on Q, there exists a Markov process on Q corresponding to (P(t)) t >o with 
initial distribution \i and cadlag sample paths. We also show that M ZtZ i jWjW i is the 
unique invariant measure for this Markov process. 

1.5. Markov process on Gelfand-Tsetlin schemes. Via the correspondence be- 
tween the probability measures on Q and central measures on paths in GT, the 
semigroup (P(t)) t >o defines a Markov evolution of central measures. It is natural to 
ask if there exists a Markov process on all probability measures on paths in GT that 
agrees with the one we have when restricted to the central measures. We construct 
one such process; let us describe its transition rates. 

Let {A^} be a starting Gelfand-Tsetlin scheme. Then 
• Each coordinate tries to jump to the right by 1 with rate 

(Af-i-* + l)(Af -i-z' + l) 
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and to the left by 1 with rate 

(X^ + k-i + w)(X^ + k-i + w'), 
independently of other coordinates. 

• If the -clock of the right jump rings but A- fc ^ = X*^ 1 , the jump is blocked. If 

its left clock rings but A,- = X^ 1 ^ , the jump is also blocked. (If any of the two jumps 
were allowed then the resulting set of coordinates would not have corresponded to 
a path in GT.) 

• If the right A^-clock rings and there is no blocking, we find the greatest number 

I > k such that = A^ for j = k, k + 1, . . . , I, and move all the coordinates 

{^i Vj=k ^° ^h e r ight by one. Given the change X\ i— > X\ k ^ + 1, this is the minimal 
modification of the initial Gelfand-Tsetlin scheme that preserves interlacing. 

• If the left Af } -clock rings and there is no blocking, we find the greatest number 
I > k such that X^-_ k = K f° r J — k,k + 1, . . . ,1, and move all the coordinates 

{A^ _ fc H =jfc to the left by one. Again, given the change A^ i-» — 1, this is the 
minimal modification of the set of coordinates that preserves interlacing. 

Since the update rule for each coordinate A- fc ^ typically depends only on a few 
surrounding coordinates, one can argue that we have a model of local random growth. 
It should be compared to the models treated in |BF08+ |, [BKlOj . where a similar 
block-push mechanism was considered with constant jumps rates, and in | BorlO+ |, 
where the jump rates were also dependent on the location and numbering of the 
coordinates. 

The key new feature of the Markov process above is the absence of the limit 
shape phenomenon. Often taken for granted in local growth models, it is simply 
nonexistent here. 

This fact becomes more apparent if we restrict ourselves to coordinates {X[ }j>i 
only. The evolution of this set of coordinates is also Markov, and it represents 
a kind of an exclusion process. Our results imply that this process has a unique 
equilibrium measure. Moreover, with respect to this measure, the asymptotic density 
linij-i.oo Xi /j is well-defined and random. It changes over time, and its distribution 
is given by a solution to the classical Painleve VI (second order nonlinear) differential 
equation, cf. |BD02| . 

1.6. Analytic continuation viewpoint. We have so far required the parameters 
(z, z', w, w') to satisfy (11.21) and (II. 3p . However, all the results would hold if (11.21) is 
replaced by more general conditions, see [ Ols03] for a precise description, with the 
only difference being that the state spaces for our Markov processes would become 
smaller. In particular, if we choose 

z = k G Z> , z' — k + a — 1, w = I e Z> , w' = I + b — 1, a, b > 0, 
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then we have to restrict ourselves to Gelfand-Tsetlin schemes with —I < < k for 
all i,j > 1. As the result, there are only k + I nontrivial parameters remaining on 
the boundary, and (P(t)) t > turns into a finite-dimensional diffusion with an explicit 
second order differential operator as its generator. The equilibrium distribution (i.e. 
the spectral zu>-measure) becomes the (k + Z)-point Jacobi orthogonal polynomial 
ensemble. See Subsection 110.11 for details. 

One can thus think of our construction as of an analytic continuation in param- 
eters (k, I, a, b) of a very well understood finite-dimensional diffusion. This point of 
view can be very fruitful: In | Olsl0+| it was heavily exploited in the construction 
and analysis of the Markov process preserving the spectral ^-measure arising from 
representation theory of the infinite symmetric group. In that case, the starting 
point for analytic continuation was the Laguerre orthogonal polynomial ensemble 
and the corresponding diffusion, rather than the Jacobi one that we have here. 



1.7. Pregenerator. As our construction of the semigroup (P(t)) t >o is fairly inex- 
plicit, it is tempting to look for its alternative definition, for example, via a generator. 

We were able to find a countable set of 'coordinates' on Q such that the action of 
the generator of (P(t))t>o on polynomials in these coordinates is given by an explicit 
formal second order differential operator, see Subsection I10.2l below. However, it re- 
mains a challenge for us to derive properties of our Markov process (or its existence) 
from the resulting formula. 



1.8. Further questions. As explained in [Ols03j, the slices GTjy can be embedded 
into the boundary so that as N — > oo, their images form an increasingly fine grid in 
Q. It is known that the A^th level z w- measures weakly converge to the spectral ones 
under these embeddings. It would be desirable to prove a similar statement for the 
Markov semigroups. 

Verifying semigroup convergence would pave the way to proving that the equilib- 
rium Markov process on the boundary can also be described as a time-dependent 
determinantal point process. The fact that the dynamical correlation functions are 
determinantal on each GUV easily follows from known techniques, although deriv- 
ing useful formulas for the correlation kernel is a separate task. Another possible 
corollary of the semigroup convergence would be that the spectral ziy-measure is a 
symmetrizing (not just an invariant) measure for (P(t)) t > . 

It seems important to continue the study of the pregenerator started in Subsection 
110.21 For example, it would be nice to understand if the space of polynomials in our 
coordinates is a core for the generator of the Markov process, and if not then how 
that space should be modified. 

Another way to benefit from investigating the generator would be to obtain a 
proof of the continuity of trajectories for our processes; we are only able to show 
that the process has cadlag trajectories at the moment. 
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All these questions and more have been settled in the case of the ^-measures 
treated in | Olsl0+ |. Unfortunately, key features of that model are not present 
here (like decomposition of the process into a one-dimensional one and a process 
on a compact set, or the existence of a convenient set of functions on the state 
space isomorphic to the well-studied algebra of symmetric functions), and one would 
clearly need new ideas. 

1.9. Organization of the paper. In Section [2] we present an abstract scheme of 
constructing a Markov semigroup on the boundary out of a consistent family of 
semigroups on the slices. In Section [3] we describe how the Gelfand-Tsetlin graph 
fits into this abstract scheme. Section H] is a brief collection of general facts about 
continuous time Markov chains on countable spaces. Section [5] provides the con- 
struction of the Markov chains on GIV's. In Section [6] we verify the consistency of 
these Markov chains. Section [7] contains a brief description of the ziy-measures. In 
Section [H] we develop a general formalism of building continuous time Markov chains 
on paths out of a consistent family of those on the slices. In Section [9] we apply this 
formalism to our specific example and discuss the exclusion type processes. Section 
[TU1 is an appendix without proofs; it contains a description of the finite-dimensional 
case of integral parameters z and w, and an explicit formula for the generator of our 
Markov process in certain coordinates. 

1.10. Acknowledgements. A. B. was partially supported by NSF grants DMS- 
0707163 and DMS-1006991. G. O. was supported by the RFBR grant 08-01-00110, 
the RFBR-CNRS grant 10-01-93114, and the project SFB 701 of Bielefeld University. 

2. Abstract construction 



2.1. Markov kernels. For a more detailed exposition, see e.g. |Mey66 Ch. IX]. 



Let E and E' be measurable spaces. A Markov kernel K : £7 — >■ E' is a function 
K(x,A), where x G E and A C E' is a measurable subset, such that K(x, • ) is a 
probability measure on E' and K( ■ , A) is a measurable function on E. 

Let B(E) and B(E') denote the Banach spaces of real- valued bounded measurable 
functions with the sup-norm on E and E', respectively. A Markov kernel K : E — > E' 
induces a linear operator B(E') — > B(E) of norm 1 via (Kf)(x) = f E , K(x, dy)f(y). 

For two Markov kernels K\ : E — > E' and K 2 : E' — > E", their composition 
K\ o K2 : E — > E" is also a Markov kernel. 

Denote by Ai ( • ) the Banach space of signed measures of bounded variation with 
the norm given by the total variation. Let be the cone of finite positive 

measures, and let Ai p { ■ ) be the simplex of the probability measures. 

A Markov kernel K : E — >• E' also induces a linear operator Ai(E) — > Ai(E') of 
norm 1 via (fiK)(dy) = J E p,(dx) K (x , dy) . This operator maps M. + (E) to M. + (E') 
and Aip(E) to Ai p (E'). Note that 5 X K = K(x, •), where 5 X is the Dirac delta- 
measure at x G E. 
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The space JA(E) (and hence Ai + (E) and M. P {E)) is equipped with a a-algebra 
of measurable sets: Any preimage of a Borel set under the map /x t— > fi(A) from 
Ai(E) to R for any measurable A is measurable. 

2.2. Feller kernels. Let P and E' be locally compact topological spaces with count- 
able bases. Let us take Borel a-algebra as the a-algebra of measurable sets for both 
of them. 

Let C( • ) C B( • ) be the Banach space of bounded continuous functions, and let 
Cq( - ) C C( • ) be its subspace of functions that tend to at infinity. 

Definition 2.1. A Markov kernel K : E — > E' is called Feller if the induced map 
B(E') ->■ B(E) maps C (P') to C (E). 

Note that different authors may use different (nonequivalent) definitions for the 
Feller property. 

The convenience of the space C ( • ) is based on the fact that this space is separable 
(as opposed to C(-) which is not separable, except in the case when the initial 
topological space is compact), and M.{ ■ ) is its Banach dual. 

2.3. Feller semigroups. A Markov semigroup is a family of Markov kernels P(t) : 
E ->• E, where t > 0, P(0) = 1 (in the obvious sense), and P(s)P(t) = P(s + t). 
Such a semigroup induces a semigroup of linear operators in B(E) as well as a 
semigroup of linear operators in Ai p (E), see above. 

We say that a Markov semigroup (P(t)) t >o is Feller if 

• E is a locally compact topological space with countable base; 

• the corresponding operator semigroup in B(E) preserves Cq(E); 

• the function t h-> P(t) is strongly continuous, i.e. t h-> P(t)f is a continuous map 
from [0, +oo) to C (E) for any / G C (E) (an equivalent condition is the continuity 
at t = 0). 

2.4. Feller semigroups and Markov processes. For more details, see e.g. [EK86, 
IV.2.7]. 

Let E be a locally compact separable metric space, and let (P(t)) t > be a Feller 
semigroup on E. Then for each /x G A'lp(-E'), there exists a Markov process corre- 
sponding to (P(t))i> with initial distribution /i and cadlag sample paths. Moreover, 
this process is strongly Markov with respect to the right-continuous version of its 
natural filtration. 

2.5. Boundary. Let E±, E2, ■ ■ ■ be a sequence of measurable spaces linked by Markov 
kernels 

A^ +1 : E N+1 -> E N , iV = l,2,... . 
Assume that we have another measurable space E^ and Markov kernels 

: — > E N , N = 1, 2, . . . , 
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such that the natural commutativity relations hold: 

A» +1 oA^ = A» N=l,2,.... (2.1) 

The kernels A^ +1 induce the chain of maps, cf. 12.11 

> M P {E N+1 ) ->• M P (E N ) -> > M P {E 2 ) -> Mp^), (2.2) 

and we can define the projective limit h_mA^ p (-EAr) with respect to these maps. By 
definition, it consists of sequences of measures (^n)n>i, £ M-p{En)i that are 
linked by the maps from ( 12. 2p . The space \imA4 p (E N ) is measurable; the cr-algebra 
of measurable sets is generated by the cylinder sets in which /i^v must lie inside a 
measurable subset of M. p (En), and all other coordinates (fik)k^N, are unrestricted. 

Observe that to any fi^ G J^i. p {E 00 ) one can assign an element of hm A4 p (E N ) by 
setting fijsr equal to the image of [i^ under the map .Mp(-Eoo) — > M. p (En) induced by 
A~. The commutativity relations ( 12. ip ensure that the resulting sequence (//jv)jv>i 
is consistent with ( 12. 2p . 

Definition 2.2. We say that E^ is a boundary of the sequence (En)n>i if the map 
A4 p (£'oo) — ► hjn A4 P (E N ) described in the previous paragraph is a bijection and also 
an isomorphism of measurable spaces. 

2.6. Feller boundary. In the setting of the previous subsection, let us further 
assume that (En) n>i and Eoo are locally compact topological spaces with countable 
bases, and all the links (A^ +1 ) 7V>1 , (A^) N>1 are Feller kernels, cf. Subsection 12.21 
Then if E^ satisfies Definition I2.2[ we shall call it the Feller boundary for (En)n>i- 

According to Subsection [221 the links (A^) N>1 induce linear operators C (E N ) — > 
Co(E OQ ). 

Lemma 2.3. The union of images of these maps over all N > 1 is dense in the 
Banach space Cq(E oc ). 

Proof. Since J^A(E 00 ) is the Banach dual to C (E oo ), it suffices to verify that if 
H G A^(-Eqo) kills all functions in our union then fi = 0. 

Assume fi is a signed measure on E^ that kills the image of Cq(En), N > 1. 
This is equivalent to saying that fiK^ = for all A^ > 1. We can represent /i as 
difference of finite positive measures 

M = an' — //, n" G MpiEvo), a,/3>0. 

Hence, a/j,'K% = /3/i"K^ for all N > 1. Since n'Kfi and fi"Kfi are in M P (E N ), 
we must have a = (3, and fi'K^ = ^"K^ . Definition 12.21 implies // = /i", thus 
fi = 0. □ 

2.7. Extension of semigroups to the boundary. In the setting of Subsection 
12. 5[ assume that for any AT > 1, we have a Markov semigroup (-P/v(£))t>o on En, 
and these semigroups are compatible with the links: 

P N+1 (t)oA% +1 = A% +1 oP N (t), t>0, N = 1,2,.... (2.3) 
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Proposition 2.4. In the above assumptions, there exists a unique Markov semigroup 
P(t) on such that 

P(t)oA% = A%oP N (t), t>0, N = l,2,.... (2.4) 

If Eoa is Feller (cf. Subsection \2.6\) and (-Pw(i))t>o a Feller semigroup for any 
N > 1, then (P(t)) t >o is also a Feller semigroup. 

Proof. Denote by 5 X the delta- measure at a point x G E^. To construct the semi- 
group (P(t))t>o, we need to define, for any t > 0, a probability measure P(t;x, ■ ) 
on Ex,. This measure has to satisfy 

P(t; x, ■ )A% = 8 X {A% o p N (t)), N>1. 

The right-hand side defines a sequence of probability measures on Erfs, and (12. ip . 
(12. 3 p immediately imply that these measures are compatible with maps (I2.2p . Hence, 
we obtain an element of hjn Ai p (E^), which defines, by definition of the boundary, 
a probability measure on E^. The dependence of this measure on x is measurable 
since this is true for any of its coordinates. 

Thus, we have obtained a Markov kernel P(t) which satisfies 

5 x (P(t) o A°°) = 5 X (A% o p N (t)), N>1, 

which is equivalent to (12 ,4p . 

To verify the semigroup property (Chapman-Kolmogorov equation) for (P{t))t>o 
it suffices to check that 

(P(s) o P(t)) oA™ = P(s + t)o A™, s,t>0, N>1, 

and this immediately follows from (12.41) and the corresponding relation for (.P/v(i))t>o- 
The uniqueness is obvious since P(t) is uniquely determined by (P(t) o A^)n>i 
that are given (12 ,4p . 

Finally, let us prove the Feller property assuming that the boundary is Feller and 
all (-P/v(t))t>o are Feller. 

We need to show that for / e Cq(E 0O ) we have P(t)f G Cq{E 00 ), and that P(t)f 
is continuous in t in the topology of Cq(Eoo). Both properties can be verified on 
a dense subset. Lemma 12.31 then shows that it suffices to consider / of the form 
/ = A™f N with f N G C (E N ). By (J22D 

P(t)f = P(t)(A%f N ) = A%(P N (t)f N ), 

which is in C (E oo ) because and -Pjv(t) are Feller. The continuity in t is obvious 
as P/v(*)/at is continuous in t, and : C (E N ) — > C (-Eoo) is a contraction. □ 

It is worth noting that our definition of the semigroup P(t) is nonconstructive: 
We are not able to describe P(t; x, A) explicitly, and we have to appeal to the 
isomorphism in Definition 12.21 instead. Thus, the difficulty in making P(t) explicit 
is hidden in the implicit nature of that isomorphism. 
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2.8. Invariant measures. In the setting of Subsection 12.51 assume that for any 
N > 1, there exists /x/v £ -M P (E N ) such that ^^P^it) = fijy (i.e., /xjv is an invariant 
measure for (-Fjv(*))t>o)- If we assume that /iAr's are compatible with the links, 

fi N+ iA% +1 = mat, iV > 1, 



then, via Definition 12.21 they yield a measure /x G -M p (Poo) such that /xA^ 1 = /xjv 
for any N > 1. Note that /x is uniquely determined by its coordinates. 
One easily sees that /x is invariant with respect to (P(£))t>o. Indeed, 

(jiP(t))A% = (jiA%)P N (t) = fi N P N (t) = fM N = fiA%. 

Moreover, if /x n is a unique invariant measure for (-P/v(t))t>o f° r an Y AT > 1 then 
the invariant measure for (P(i)) t > is unique too as its convolution with must 
coincide with /xjv- 

3. Specialization. Gelfand-Tsetlin graph 

3.1. Spaces and links. Let A^ be a positive integer. A signature A of length N 
is an A^-tuple of weakly decreasing integers: A = (Ai > • • • > A a?) G Z N . Denote 
by GTat the set of all signatures of length A" (the notation GT is explained below). 
This countable set will serve as our space E N from the previous section. 

Signatures of length A" parameterize irreducible representations of the unitary 
group U(N) and are often referred to as highest weights, cf. |Wey39| , |Zhe70] . For 
A G GTat denote the corresponding representation by 7i\, and denote by DiniAr A the 
dimension of the corresponding linear space. It is well known that 

Ili<i<j<N( X i ~ i ~ X J + h 



Dim N \ = — jrr-- : , A G GT N . 



Define a matrix [A^ +1 (A, v)\ A6GTjv+l ueG j N with rows parameterized by GTtv+i 
and columns parameterized by GTat via 

, IIi<i<j<jv(^-^-^-+j) A 

as +i (a,i/) = <{ ''n^w^-^-^+i)' 1 

0, otherwise, 
where the notation v -< A stands for interlacing: 

^ -< A •<=>- Ai > i>i > A 2 > ^2 > ■ ■ ■ > ^ > Ajy +1 . 

, \~^- /-trim r\ 1 ct/~\ V\ /-» urn -^--i - i^m l vi -i - In , 



Note that the nonzero entries of A^ +1 can also be written in the form 



«') = S. (3-D 

It is not hard to show that Aj^ +1 is a stochastic matrix: ^j/ggtjv -^-iv +1 (^' ^) = 1 
for any A G GTjv+i- Indeed, DmiAr+x A is equal to the number of the sequences 
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(known as Gelfand-Tsetlin schemes, thus the notation GT) 

A (D _< a( 2 ) -< • ■ • -< X^ = X, A^ e GTj, 

and Aj^ +1 (A , v) is the fraction of the sequences with X^ = v. The stochasticity 
also follows from the branching rule for the representations of unitary groups: For 
any A G GT/v+i, 

K\\u(N) ~ (J) Tfv 

The matrices A^ +1 viewed as Markov kernels A^ +1 : GIV+i — > GT^r are our 
links, cf. Subsection 12.51 Set GT = |J Ar>1 GTAr. We endow GT with the structure 
of a graph: Two vertices A and v are joined by an edge if and only if v -< A or A -< v. 
This graph is called the Gelfand-Tsetlin graph, and the matrix elements of the links 
are often called cotransition probabilities for this graph, cf. |Ker03] . 

3.2. Boundary, let IR+ C K be the set of nonnegative real numbers and be the 
product of countably many copies of M.+. Consider the space 

R 4oo+2 . = ^oo x R oo x R oo x R oo x R+ x R+ 

and equip it with the product topology. We choose to be the closed subset 
Q C M^°° +2 consisting of the sextuples 

u = (a + , (3 + , a~, (3~ , 5 + , <T) e M.^° +2 

satisfying the conditions 

a ± = (af>af >...), /3± = (fif > fit > ■ ■ ■ ), ^ > 0, 

oo 

^K ± + /3f)<5 ± , # + 
i=i 

One easily sees that Q is a locally compact metrizable topological space with a 
countable base. We endow Q with the corresponding Borel structure which makes 
Q a measurable space. 

It will be convenient to use the notation 

oo 

i=i 

Define the projections/links : Q — > GTjv, N > 1, by 

A^(w,A) = Dim JV A-det[^ Ai _ i+i ]y =1 , w G fi, A e GT N , (3.2) 
where {<f n }n=-oo are the Laurent coefficients of the function (\u\ = 1) 

$ ._ r 7+(^D+7-(n- 1 -l) TT 1 + Pt( U - 1) 1 + ~ 1) _ n 

*" ( " ) - fi l - «, + (« - i) i - «■(«- - i) " „^/"" ' 

(3.3) 
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Theorem 3.1. The space Eoo = Q is the boundary of the chain of spaces (En = 
GT N ) N >i with links as above in the sense of Definition \2. 6 A 

Proof. This result is essentially proved in §9 of |Ols03j ; we provide below some 
necessary additional comments. 

Let us abbreviate A N = A4 P (GT N ) and A = hm A4 P (GT N ) = \^mA N . As in 
|Ols03] we embed A into the vector space of all real- valued functions on the set of 
vertices of the Gelfand-Tsetlin graph. That space is endowed with the topology of 
pointwise convergence, and A inherits this topology. The measurable structure of 
A is the Borel structure corresponding to this topology. 

Obviously, A is a convex set; let Ex A C A denote the subset of extreme points. 
Theorem 9.2 in [Ols03j (which is based on Choquet's theorem) says that Ex A is a 
Borel subset of A, and each point of A is uniquely representable by a probability 
Borel measure concentrated on Ex A. On the other hand, it is readily seen that, 
conversely, any probability Borel measure on A (in particular, on Ex A) represents 
a point of A, the barycenter of that measure. This gives us a bijection between 
.Mp(ExA) and A. 

The next step consists in identifying the abstract set Ex A with the concrete 
space Q. This is achieved with the help of Theorem 1.3 in |Ols03] . Namely, as 
is pointed out in the proof of Theorem 9.1 in [Ols03j . there is a natural one-to- 
one correspondence between the points of Ex A and the extreme characters of the 
infinite-dimensional unitary group U(oo), which in turn are parameterized by the 
points of the space Q, see Theorem 1.3 in |Ols03] . 

Then we have to verify that the embedding Q — > A induced by the identification 
Q = Ex A is given by the kernels A^. This is shown by the computation in |Voi76] . 

We have thus constructed a bijective map A4 P (Q) — » A, and it remains to prove 
that it is a Borel isomorphism. As shown in the proof of Theorem 8.1 of |Ols03] . 
the map oj h-» A^(cj, A) is continuous for every N = 1,2,... and every A G GT N . 
This implies that the map Ai p (Q) — > A is Borel. To show that the inverse map 
is also Borel one can apply an abstract result (Theorem 3.2 in |Mack57j ). which 
asserts that a Borel one-to-one map of a standard Borel space onto a subset of a 
countably generated Borel space is a Borel isomorphism. This result is applicable in 
our situation, since the Borel structure of Q is standard, so that the induced Borel 
structure on Ai p (Q) is standard, too. 

□ 

Remark 3.2. Observe that the maps on (GTat)tv>i consisting of shifts of all coor- 
dinates of signatures by 1, 

A = (Ai, . . . , Ajv) A = (Ai = Ai + 1, . . . , X N = X N + 1), 

leave the links intact: A^ +1 (A,i/) = A^ +1 (A, v). There is also a corresponding 
homeomorphism of Q, which amounts to the multiplication of the function $ w (w) 
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by u: For u = (a ± , fi±, 5 ± ) G fi define £5 = (5 ± , /3 ± , 5 ± ) G fi by 

fit = i- ^, (Pt, ft,---) = (/3+,/3+,...), (/3r,&,...) = (^>&v-0 

(note that > fit because fit + fii < 1)- Then (13. 2 p and the relation 
u(l + ^-(u- 1 - 1)) = 1 + (1 - fi{)(u - 1) 

show that Afi{u), A) = A~(u3, A) for any A G GTjy and N > 1. 

This automorphism of the Gelfand-Tsetlin graph and its boundary has a repre- 
sentation theoretic origin, cf. Remark 1.5 in |Ols03] and Remark 3.7 in [BO05aJ. 



3.3. The boundary is Feller. Following definitions of Subsection 12.61 in order to 
show that = Q is a Feller boundary of the chain (E N = GTtv)at>i we need to 
verify two statements: 

• the spaces (En)n>i and are locally compact topological spaces with count- 
able bases; 

• the links (A^ +1 )at>i and (A^)at>i are Feller kernels. 

The first statement is obvious from the definitions. The goal of this subsection is 
to prove the second one. 



Proposition 3.3. For any N > 1, the linear operator i3(GlV) — > i3(GlV+i) in- 

V+i)- 



duced by the Markov kernel A^ +1 maps Cq(GTn) to Co(GT/ 



Proof. As the norm of the linear operator in question is equal to 1 and Co ( • ) is a 
closed subspace of B( • ), it suffices to check that the images of all delta-functions on 
G1V are in Co(GTat + i). 

For a v G GT^, let 8 V be the delta-function on GTat concentrated at v. Then for 
A G G1V +1 

Ul<i<j<N(^ ~ 1 - + J) ., jX 

if v -< A, 



{A N N +l 5 u ){\) = { ■ rii^W^ - i - A; + j) ' 

otherwise. 

If we assume that (A^ +1 5 U )(X) ^ then A — > oo is equivalent to either Ai — >■ +oo, 
or \n+i — > — oo, or both; all other coordinates must remain bounded because of the 
interlacing condition v -< A. But then it is immediate that at least one of the factors 
in the denominator in (A^ +1 6 u )(\) tends to infinity. Thus, for any fixed v G GUV, 
(A^ +1 S V )(X) ->■ as A -> oo as needed. □ 

Proposition 3.4. For any N > 1, the linear operator i3(GTjy) — > B(Q) induced by 
the Markov kernel A^ maps Cq(GTjv) to Cq(Q). 



Proof. As in the proof of Proposition 13.31 it suffices to prove that for any v G GTjv, 
(Atf5 u )(w) = Ax(u),\) belongs to Cq(Q) as a function in oj G Q. The proof of 
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continuity of A^(uj, A) in u is contained in the proof of Theorem 8.1 of [Ols03j. It 
remains to show that A^(oj, A) — > as uj — > oo. Note that uj — >• oo is equivalent to 
5 + + 5~ +oo. 

Observe that the coefficients ip n = (p n (uj) from (13.31) can be written as 

I f du 

<p n (u) = — i QM-^z v neZ - ( 3 - 4 ) 

2m J H=l u n+l 

Note that < 1 on the unit circle \u\ = 1, because the modulus of each of the 

factors (1 - af{u ±l - l))" 1 , (1 + Pf(u ±l - 1)), and e^^" 1 ) is < 1. 

We are going to prove that for any fixed n G Z, ip n — > as uj — > oo; by (13. 2 p this 
would imply the needed claim. 

Let us assume the converse, i.e. assume that there exist c > 0, n G Z, and a 
sequence {a;(fc) } fc > x C Q with limfc_ >tX) a;(A;) = oo, such that ip no (cj(k)) > c. Let 
us denote by (a ± (fc), /3 ± (/c), <5 ± (/c)) the coordinates of uj(k). We will now collect 
information about {uj(k)} that will eventually lead to a contradiction. 

Step 1. We must have sup fc>1 af (k) < oo. Indeed, if there is a subsequence {k m } m >i 
such that af(k m ) — > oo, then along this subsequence (1 — af(u ±l — tends to 
zero uniformly on any compact subset of {u : \u\ — 1} \ {u — 1}, which implies that 
the right-hand side of (13. 4p tends to zero. 
Let us fix A > such that sup fc af(k) < A. 

Step 2. Assume to ranges over the subset of elements of VL with af < A and /3f < |. 
Then for any e > 0, 

lim &uj(u) = uniformly on {u G C : \u\ — 1, Km < 1 — e}. 
Indeed, for u on the unit circle with Km < 1 — e we have elementary estimates 

|1 + p(u - 1)| 2 = (1 - /3) 2 + /3 2 + 2/3(1 - /3)Kw 

= 1 - 2/3(1 - /3)(1 - Km) < 1 - 2/3(1 - /3)e < 1 - 0e < e~ Pe , (3.5) 

|1 - a{u - 1)T 2 = (1 + 2a(l + a)(l - Km)) -1 

< (1 + 2a(l + a)e)- x < (1 + 2ae)- 1 < e - constae , (3.6) 

, 7 +( u _l) +7 -( u -l_l)| 2 _ e _2( 7 ++ 7 -)(l-SR«) < e -2(7++7-)e^ 

with a suitable constant const > (that depends on A). Thus, if 

oo 

6 + + <T = 7 + + 7" + + $ + + &~) ~+ 00 

i=i 
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then at least one of the right-hand sides in these estimates yields an infinitesimally 
small contribution, and ^^(u) must be small. Thus, under the above assumptions 
on oj, we see that u — > oo implies (p n (u) — > uniformly on n e Z. 

Step 3. Now we get rid of the restriction f3f < |. Set 

B ± (k) = #{* > 1 | Pt(k) > \}. 

Since for any k > 1 we have f3f{k) + Pi{k) < 1, at least one of the numbers B ± (k) 
is equal to 0. The statement of Step 2 shows that for any subsequence {oJk m } of our 
sequence {u(k)}, we must have B + {k m ) +B~(k m ) — > oo. Hence, possibly passing to 
a subsequence and switching + and — , we may assume that B + (k) — > oo as k —¥ oo. 
Using the identity (cf. Remark 13. 2p 

1 + /3{u - 1) = w(l + (1 - Z?)^- 1 - 1)) 

B + {k) times on ^(^(it), we see that y no (a;(A;)) = y> no (a) (&)), where Cj{k) 
is obtained from o;(fc) as follows: Each /3 + -coordinate of u(k) that is > 1/2 is 
transformed into a (3~ coordinate of u(k) equal to 1 minus the original /3 + -coordinate; 
all other coordinates are the same (equivalently, the function $o>(ifc)("u) is multiplied 
by u~ B( V). Let (^(k), P ± (k),^ ± (k) 1 5 ± (k)) be the coordinates of u(k). 

Step 4- Since no /3-coordinates of u(k) are greater than 1/2, the argument of Step 
2 implies that if 5 + (k) + 5~(k) — > oo then ip m (u(k)) = ip no -B+(k){&(k)) — >• as 
k — > oo, which contradicts our assumption. Hence, 5 + (k) + 5~(k) is bounded. 

Let us deform the integration contour in f !3.4p to \u\ — R with A/(1+A)<R <1. 
Using the estimates (for |u| = R, < a < A, < < 1/2) 

|1 + ${u ±l - 1)| < 1 + p\u ±l - 1| < e constl/3 , 

|1 - aiu^ 1 - I)]' 1 < |1 - a^ 1 - 1)T X < e^ 11 ^", 

| e 7(« ±:L -l)| ^ e const 3 7 

with suitable constj > 0, j = 1,2,3, we see that |$o( fc )(u)| < e const4 ^ +(fe)+ ^ (fc)) , for 
a const4 > 0, which remains bounded. On the other hand, the factor M - n o-i+-B+(fc) 
in the integral representation H3.4[) for (p no _ B +^(u;(k)) tends to uniformly in u, 
\u\ = R < 1. Hence, ip no (uj(k)) = (p m -B+(k)(ti(k)) — > as k — > oo, and the proof of 
Proposition 13.41 is complete. □ 

4. Generalities on Markov chains on countable spaces 

4.1. Regularity. Let E be a countable set, and let (P(t)) t >o be a Markov semigroup 
on E 1 . Each P(t) may be viewed as a matrix with rows and columns marked by 
elements of E; its entries will be denoted by P(t;a,b), a,b e E. By definition, 
P(t; a, 6) is the probability that the process will be in the state b at the time moment 
t conditioned that it is in the state a at time 0. Thus, all matrix elements of 
P(t) are nonnegative, and their sum is equal to 1 along any row - the matrix P(t) 
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is stochastic. The transition matrices P(t) also satisfy the Chapman-Kolmogorov 
equation P(s)P(t) = P{s + t). 

Assume that there exists an E x E matrix Q such that 

P(t;a,b) = l a=b + Q(a,b)t + o(t), t±0. (4.1) 

This relation implies that Q(a, b) > for a ^ b and Q(a, a) < 0. Further, we will 
always assume that 

Q(a, b) = —Q(a, a) for any a G E. 

b+a 

This is the infinitesimal analog of the condition YlbeE ci,b) = 1. 

It is well known that the Chapman-Kolmogorov equation implies that P(t) satis- 
fies Kolmogorov's backward equation 

±P(t) = QP(t), t>0, (4.2) 

with the initial condition 

P(0)=ld. (4.3) 

Under certain additional conditions, P(t) will also satisfy Kolmogorov's forward 
equation 

j t P(t) = P(t)Q, t>0. (4.4) 

One says that Q is the matrix of transition rates for (P(t)) t > . 

One often wants to define a Markov semigroup by giving the transition rates. 
However, it may happen that this does not specify the semigroup uniquely (then 
the backward equation has many solutions). Uniqueness always holds if E is finite 
or, more generally, if E is infinite but the diagonal entries Q(a, a) are bounded. 
However, these simple conditions do not suit our purposes, and we need to go a 
little deeper into the general theory. 

Let us write Q in the form Q = —q + Q, where — q is the diagonal part of Q and 
Q is the off-diagonal part of Q. In other words, 

q(a,b) = -Q(a,a)l ab , Q(a,b) = < 

10, a = b. 

Define P^ n \t) recursively by 

pl°](t) = e~ tq , P [n] {t) = [ e~ Tq QP [n ' 1] (t -r)dr, n > 1, 

Jo 

and set 

oo 

P(t) = ^pN( t ) ) t > . 

n=0 
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Theorem 4.1 (|Fel40j). (i) The matrix P(t) is substochastic (i.e., its elements are 
nonnegative and ^ b P(t]CL, b) < 1). Its elements are continuous in t G [0, +oo) and 
differentiable in t G (0, +oo), and it provides a solution of Kolmogorov's backward 
and forward equations (" 14 .21) . ( I4.4p with the initial condition ( 14.31) . 
(ii) P(t) also satisfies the Chapman- Kolmogorov equation. 

(Hi) P(t) is the minimal solution of ( 14. 2 p (or ( 14. 4p ) in the sense that for any 
other solution P(t) of ( 14. 2 p (or ( 14 - 4 j) ) with the initial condition ( I4.3P in the class of 
substochastic matrices, one has P(t; a, b) > P(t; a, b) for any a,b G E. 

Corollary 4.2. If the minimal solution P(t) is stochastic (the sums of matrix ele- 
ments along the rows are all equal to 1) then it is the unique solution of (14. 2 p (or 
( 14. 4p ) with the initial condition (14. 3 p in the class of substochastic matrices. 

If the minimal solution P(t) is stochastic one says that the matrix of transition 
rates Q is regular, cf. Proposition 14.31 

Observe that the construction of P(t) is very natural: the summands P^(t; a, b) 
are the probabilities to go from a to b in n jumps. The condition of P(t) being 
stochastic exactly means that we cannot make infinitely many jumps in a finite 
amount of time. 

A much more detailed account of Markov chains on countable sets can be found 
e.g. in |And91] . 

Later on we will need the following sufficient condition for P(t) to be stochastic. 

For any finite X, X C E, a G X, denote by T at x the time of the first exit from 
X under the condition that the process is in a at time 0. Formally, we can modify 
E and Q by contracting all the states b G E \ X into one absorbing state b with 
Ql c = for any c G X U {b}. We obtain a process with a finite number of states for 

which the solution P(t) of the backward equation is unique. Then T a x is a random 
variable with values in (0, +oo] defined by 

Prob{T a , x <t} = P(t;a,b). 

Proposition 4.3. Assume that for any a G E and any t > 0, e > 0, there exists 
a finite set X(e) C E such that Prob{T a x(e) < t} < e. Then the minimal solution 
P(t) provided by Theorem \4-l\ is stochastic. 



Proof. Consider the modified process on the finite state space X(e) U {b} described 
above. Since its transition matrix P(t) is stochastic, 

P(t] a,b) = l- P(t- a,b) >l-e. 

The construction of the minimal solution as the sum of P'^'s, see above, immediately 
implies that P(t; a, b) > P(t; a, b). Thus, Y2b -^"(^ a,b) > 1 — e for any e > 0. □ 
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4.2. Collapsibility. In what follows we will also need a result on collapsibility or 
lumpability of Markov chains on discrete spaces. Let us describe it. 

Let E = \_\ ieI Ei be a partition of the countable set E on disjoint subsets. Assume 
we are given a matrix Q E of transition rates on E and a matrix Qj of transition 
rates on / such that 

^ Q E (a,b) = Qi(i,j) for any aeE h i,jeI. (4.5) 

beEj 

Denote 

q E (a) = -Q E (a,a), qj{i) = -Qi(i,i). 
For any i E I, let Qj be a matrix of transition rates on Ei defined by 

Qi(a, b) = Q E (a, b) if a ^ b, a,b G E u 
qi(a) = -Qi(a,a) = ^ Qi(a,b), a E Ei. 

b€Ei,b^a 

Observe that q E {o) = qi{o) + qi(i) for any a G Ei. 

Denote by P E (t) and Pi(t) the minimal solutions of Kolmogorov's equations for 
Q E and Qj, respectively. 

Proposition 4.4. Assume that for any i e I, Qi is regular. Then for any t > 

J2PE(t;a,b) = Pj(t;t,j) (4.6) 

for any i,j G / and any a & Ei. In particular, if Qj is regular then so is Q E and 
vice versa. 

Proof. Let us use notations q, Q, and P^ for the diagonal and off-diagonal parts of 
the matrices of transition rates, and for the nth terms in the series representations 
of minimal solutions, respectively. 

The hypothesis means that the identity 

oo 

^^Pf ] (t;a,6) = l, aeE t , i G /, 

beEi n=0 

holds, where 

Pl n] (t) = J e- tiqi Qie-( t2 - tl)qi Qi ■ ■ ■ e-^^-^Q^-^dh ■ ■ ■ dt n 

0<ti<-<t n <t 

and Pi°\t) = e~ tq \ Using the fact that <&( • ) = q E { ■ ) — qi(i) on Ei, rewrite this 
identity as 
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oo „ 

EE / { e ~ tiqE1 E, Qe Ie^-^U ■ ■ ■ 



l Ei Qe l Ei e~ {t - tn)qE ) (a, b)dh ■ ■ ■ dt n = e~ tq ^ . (4.7) 



n=0 0<tx<-<t n <t 



The probabilistic meaning of this formula is that the time that the minimal so- 
lution -P_e(t) started at a G spends in Ei is exponentially distributed with rate 
qi(i), independent of a. 

The minimal solution Pi(t) has the form 

oo „ 

PjfaiJ) = E J E e' s ^Qj{tM)e- {s2 - Sl)qi{kl) --- 

n=0 0<si<—<s n <t fc l.---. fc n-lG-f 

• • • Qi{k n -uj)e- (f - sM) ds 1 . . . ds n . (4.8) 

By <m, 

Qi(k ) l) = Y J QE(c 1 d) for any k,l £ I, k 7^ I, c G E k . (4.9) 

d£ El 

Substituting the right-hand side of (14.91) for each Qi(- , • ) and the left-hand side of 
(14. 7p for each e~ sgi (' \ in the nth term we obtain the part of the series for P#(t; a, b) 
with a G Ei that takes into account trajectories whose projections to / make exactly 
n jumps, and in addition to that there is a summation over b G Ej. Clearly, the 
summation over n reproduces the complete series for Pe{^\ &j b) thus proving (14.61) . 

The equivalence of stochasticity of PE{t) and that of Pi(t) immediately follows 
from summation of (14.61) over j G /. □ 

4.3. Infinitesimal generator. The last part of the general theory that we need 
involves generators of Markov semigroups. 

Assume that we have a regular matrix of transition rates Q. Let (P(t)) t >o be the 
corresponding Markov semigroup and assume in addition that it is Feller. 

The generator A of the semigroup (P(t))t>o is a linear operator in Cq{E) defined 

by 

Af= i im P ®f ~ ? , (4.10) 

The set of / G Cq(E) for which this limit exists (in the norm topology of Cq(E)) is 
called the domain of the generator A and denoted by D(A). It is well known that 
the operator A with D(A) as above is closed and dissipative. 

It turns out that the domain D(A) can be characterized by an apparently weaker 
condition, which is easier to verify in practice: 
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Proposition 4.5. If f G Cq(E) is such that the limit in the right-hand side of 
( 14.1 Op exists pointwise and the limit function belongs to Cq{E), then f G D(A), so 
that the limit actually holds in the norm topology. 

Proof. The idea is that the set of couples of vectors (/,<?) G Co(E) x Cq(E), such 
that g is the pointwise limit of the right-hand side of (I4.10p . serves as the graph of 
a dissipative operator A extending A, whence A = A. A detailed argument can be 
found in |Ito06| §4.8]. In fact, [Ito06j considers the case of a compact state space 
E. However, the proof goes through word-for-word; the only property one needs is 
that for any / G Cq(E), f attains its minimum if it has negative values. □ 

The following statement is probably well known but we were not able to locate it 
in the literature. 

Proposition 4.6. Assume that for any a G E the set of b such that Q(a,b) ^ is 
finite. Then 

D(A) = {fe C (E) | Qf G C (E)}, (4.11) 
and for f G D{A), Af = Qf. 

Proof. First of all, due to the assumption on the matrix Q, Qf is well defined for 
any function / on E. We will show that for any / G Cq(E) and a G E 

lim r l Y)(P(f; a, b) - l a=b )f(b) = V Q(a, b)f(b). (4.12) 

beE beE 

Then the claim of the proposition will follow from Proposition 14.51 
Set 

X = {a} U {b' G E | Q(a,b') > 0}. 
By our hypothesis, this set is finite. We will show that 

P(t;a,b) = 0(t 2 ) } t->0. (4.13) 

b&E\X 

This would imply that we can keep only finitely many terms in (j4.12p . and then 
(KTflf would follow from (143]) . 

Observe that the left-hand side of (I4.13P is the probability of the event that the 
trajectory started at a is outside of X after time t. In order to exit X the trajectory 
started at a needs to make at least two jumps. Assume that the first two jumps 
are a — > a' — > a" with a' G X. Since X is finite, the rates of leaving a' (equal to 
—Q(a',a')) are bounded from above, and the probability of leaving X after time t 
can be estimated by 

-Q(a, a) max(-Q(a', a')) • t 2 + o{f) = 0{t 2 ) } t -> +0, 



as required. 



□ 
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Corollary 4.7. Under the hypothesis of Proposition \4~1^ assume additionally that 
for any b G E the set of a G E with Q(a, b) ^ is finite. Then any finitely supported 
function f on E belongs to D(A). 

Proof. Indeed, this follows immediately from Proposition 14. 6[ since Qf is finitely 
supported and hence belongs to Cq(E). □ 

5. Semigroups on GT n 

The goal of this section is to define Markov semigroups (-P/v(t))t>o on En = G1V 
and prove that they are Feller. 

5.1. Case N — 1. Birth and death process on Z. Let (u,u') and (v,v') be two 
pairs of complex numbers such that (u + k)(u' + k) > and (v + k)(v' + k) > 
for any k G Z. The condition on (u, u') means that either u' = u G C \ IR or there 
exists k G Z such that k < u, u' < k + 1; the condition on (v,v') is similar. Note 
that u + u' G IR and i)+fl'6l. Assume additionally that u + u' + v + v' > —1. 

Define a matrix of transition rates [T>(x, y)] x yeZ with rows and columns param- 
eterized by elements of E\ = GT\ = Z by 

(x — u)(x — u'), if y — x + 1, 

(x + u )(x + f'), if |/ = x — 1, 

— (a; — u)(x — u') — (x + v) (x + v'), if y = x, 

0, otherwise. 



V(x,y) 



(5.1) 



In the corresponding Markov chain the particle would only be allowed to jump by 
one unit at a time; such processes on Z>o are usually referred to as birth and death 
processes, while our Markov chain is an example of so-called bilateral birth and death 
processes which were also considered in the literature, see e.g. |Fel57l Section 17], 
[Pru63j . [Yan90j . 

Note that T>(x, x ± 1) > for all x G Z, because of the conditions imposed on the 
parameters. 

Theorem 5.1. The matrix of transition rates T> is regular. Moreover, the corre- 
sponding Markov semigroup is Feller. 

In what follows we denote this semigroup by (Pi(t)) 4 > . 

The proof of Theorem 15.11 is based on certain results from |Fel59] ; let us recall 
them first. 

Consider a birth and death process on Z> with transition rates given by 



Q{x,y) 



(3 X , ify = x + 1, 

S x , ify = x- 1, 

-Ac - S x , ify = x, 

0, otherwise. 
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Here {fl x }x>o, {S x }x>i are positive numbers, and we also set 5o = 0. 
The natural scale of the process is given by 

x = 0, Xk = y^ ^ A; = 1, 2, ... , = lim x k . (5.2) 

^ Po • • • Pi fc ->°° 

Note that may be infinite. Denote by A the operator on the space of functions 
on A = {xo, xi, . . . } defined by 

(Af)(xi) = -(5 t + Pi)f(xi) + 6 l f(x l . 1 ) + 0if(x i+1 ), i = 0,l,.... 

Fix n > 0. Let Fi(t) be the probability that the process started at % reaches n 
before time t. Let G n (t) be the probability that the process started at n reaches 
before time t and before the process escapes to infinity. 

Theorem 5.2 ([Fel59]). (i) For any a > there exists exactly one function u on 
A such that Au = au, u(xq) = 1. The function u is strictly increasing: u{xq) < 
u(xi) < u(x 2 ) < . . . and satisfies 

n-l 

u(x Tl ) = 1 + a ^ u{x k )(x n - x k )fi k (5.3) 



fc=0 



with 



Furthermore, 



^ = — x j—, k = l,2,..., A*o = l- (5-4) 

oi ■ ■ ■ o k 



u{Xi) 



u{x n ) 

(ii) With u(-) as above, set 



oo 

at 



e- at dFi(t), 0<i<n. (5.5) 



oo 

v{x n ) = u{x n ) > r, n = 0, 1, 

M(x J )M(x i+1 ) 



This is a strictly decreasing function, and 

= e" a ^G n (t), n = 1, 2, . . . . (5.6) 

Furthermore, \im n ^ 00 v(x n ) = if Xoo = oo and ^ n x n \i n diverges. 

The following statement is contained in Feller's paper as well, but not explicitly; 
for that reason we formulate it separately. 

Corollary 5.3. // then lim n ^ 00 -u(x n ) = oo. 
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Proof. Let us estimate the sum in the right-hand side of ( 15. 3ft : 

n—1 n—1 n—1 n 

k=0 k=0 k=0 l=k+l 

n l—l n—1 I n—1 

= 5^ y~]( x i ~ xi-dvk = ^2( x i+i - x i) 5^ ^ - _ = Xn - ( 5 - 7 ) 

1=1 k=0 1=0 k=0 1=0 

Hence, u(x n ) > 1 + ax n , and the statement follows. □ 

Let us now apply Feller's results to our situation. 

Proof of Theorem \5.1[ By Proposition I4.3[ in order to show that the minimal solu- 
tion is stochastic it suffices to prove that the probability that the first passage time 
from to n is below a fixed number, converges to zero as n — > +oo. Indeed, as 
shifts x — > x + const and sign change x i— > —x keep our class of processes intact, 
similar convergence would automatically hold for passage times to the left, and also 
for passage times from any initial position. Denote the first passage time from to 
n by T n . 

A simple coupling argument shows that T n stochastically dominates the first pas- 
sage time from to n for the birth and death process on Z> with the same transition 
rates (see (I5.ip ). except that the jump from to —1 is forbidden. Let us denote this 
new first passage time by T n . Thus, 

Prob{T n < t} < Prob{f n < t} for any n > and t > 0. 

For the application of Theorem 15.11 we then set 

Px — ( x — u) (x — v!) , x > 0, 5 X = (x + v)(x + v'), x > 1, 5 = 0. 

As 

^4 = const ^ + i + + 1 ± 1L „ cons t ■ ^'W I _ oo, 

Po---pi r{-u + l + l)T{-u' + l + l) 

our original assumption u + u' + v + v'>— 1 implies 

x k ~ const • k u+u ' +v+v ' +1 , k^oo, (5.8) 
cf. (15.21) . and Xoo = limk^ooXk = oo. Therefore, Corollary 15.31 yields 

lim u(x n ) = oo. 



On the other hand, from (I5.5P with i = and any a > we obtain 
— = e- ar dF (T)> / e~ ar qIF q {t) > e~ at qIF {t) = e~ at Prob{T„ < t] 

X n) JO JO JO 



U{X 

whence 



Prob{T n < t} < — — r ^ 0, n -> +oo. 
u{x n ) 



MARKOV PROCESSES 



27 



Since T n is dominated by T n , we have shown that our Markov chain does not run 
away to infinity in finite time, and hence it is uniquely specified by the transition 
rates. Let (Pi{t))t>o be the corresponding semigroup. 

We now need to prove that (Pi(t)) t >o is Feller. This is equivalent to showing that 
lim n ^ ±00 P\{t] n,i) — for any « G Z and t > 0. 

Shift and sign change invariance (see the beginning of the proof) imply that it 
suffices to consider i = and n — > +oo. Observe that Pi(t; n, 0) cannot be greater 
than the probability that the first passage time from n to is not more than t. Let 
us denote this first passage time by S n ; we have Pi(t; n,0) < ProbjS^ < t}. 

This first passage time is the same for our birth and death process on Z and for 
its modification on Z>o that was used in the first part of the proof. On the other 
hand, for the process on Z>o the Laplace transform of S n is given by (15. 6p . 

By (15 .4p we have, as k — > oo, 

go — gfc-l , T(-U + k)T(-u' + k) ,-2-u^-v-v' 

= -F 7— = const w , . , 1W , , , i -i \ ~ const • k 

Ox ■ ■ ■ Of. T(v + k + l)T(v' + k + 1) 

Hence, cf. (Oj) 

Xkf^k ~ const • k~ x , k — > oo, 
with a nonzero constant, and J2 n x n fi n diverges. Theorem I5.2( ii) then gives 

lim v(x n ) = 

n— >oo 

and using (I5.6P and estimating the Laplace transform as above we obtain 

Prob{£„ <t}< e °y x ") _> o, n ->• oo. 
f(x ) 

As P\{t; n, 0) < Prob{,S n < t}, the proof of Theorem 15.11 is complete. □ 

5.2. The case of general N. Let iV > 1 be a positive integer, and let (u,u') and 
(v,v r ) be as in Subsection 15.11 

Define a matrix [l>W(A, ^)] A ueGT of transition rates with rows and columns 
parameterized by points of £7/v = GUV via 

V w (X,u) = ^i)lfc=n^i> + V(l 2 , n 2 )l {h=ni ,^ 2} + ... 

+ D{l N , n N )l{ h=rHii ^ N} ) - d N ■ l x=v (5.9) 

with lj = Xj + N — j, rij — i/j + N — j, 1 < j < N, matrix T>( ■ , ■ ) as in (15.11) . and 

N(N — 1)(N — 2) , , ^N(N-l) 
d N = — - (u + u' + v + v') — — -. (5.10) 
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In other words, an off-diagonal element V^ N \X, v) can only be nonzero if there 
exists exactly one index i such that — Aj = ±1 while for all other indices j we 
have Xj = Vj. Under this condition 



V) 



(k - u)(k - u') n lj \ if^-Ai = l, 



With this explicit description, the diagonal entries of T>^ N ' have to be defined by 
V {N \X,X) = - V {N) AG GUV (5.11) 

The fact that (15.111) holds for T>^ N > defined by (15.91) will be proved in Step 1 of the 
proof of the following theorem. 

Theorem 5.4. The matrix of transition rates P'*' is regular. The corresponding 
semigroup (-Pjv(i))t>o has the form 

P N (t; A, u) = e-^ ^4tI det [ p i (*! K + N-i^ + N- j)] f A, u G GT N , 

(5-12) 

with (Pi(t))t>o as in Subsection \5.1[ Moreover, this semigroup is Feller. 
Proof. The proof of Theorem 15.41 will consist of several steps. 

Step 1. Let us show that with definition (15. 9p . relation (15.111) holds. It is convenient 
to encode signatures of length N by iV-tuples of strictly decreasing integers via 

A = (Ai > A 2 > • • • > X N ) < — > {li > l 2 > ■ ■ ■ > l N ), lj = \j+N-j, 1 < j < N. 

This establishes a bijection between GUV and the set 

%N = {{xi, • • • ) %n) ^ ^ N \ Xi > X 2 > ■ ■ ■ > X N }. 

In Xn, matrix T>^ N ' from ( 15 .9p takes the form 

rt N \X,Y) = ^|^(p(x 1) y 1 )l K=yi) ^ 1} + D(x 2 ,y 2 )l K=2/i) ^ 2} + ... 

+ V(x N ,y N )l[ Xi=yii i^ N y) -d N l x =Y (5.13) 
with X = (xt, . . . , x N ) G Xat, Y = (yi, ...,y N )e X N , and 

V N {z x , ...,z N )= ] ] (zi - Zj). 

l<i<j<N 

In this notation, (15.111) is equivalent to 

(©! + ••• + V N )V N (X) = d N V N (X), X G Z N , (5.14) 
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where T> 3 denotes a linear operator on 7* N with 

v ji x , Y ) = 

Indeed, both sides of ( 15.14ft are skew-symmetric, and restricting to Xn yields ( 15.11ft . 
Let A and V be the standard forward and backward difference operators on Z: 

A/(x) = f(x + 1) - f{x), V/(x) = fix) - fix - 1) 

for any function / : Z -»• C. Note that AV = A - V. 

One easily checks that the operator T> with matrix (15.11) has the form 

V = ctAV + tA (5.15) 

with 

a — ix + v) (x + v'), t = sx + (««' — vv'), s = — (u + u' + v + v'). 

Hence, for any m = 0, 1, 2, . . . 

T>x m = (m(m — 1) + sm) ■ x m + lower degree terms, (5.16) 

in paticular, T> preserves the degree of a polynomial. This implies that the left- 
hand side of (15.14ft is a skew-symmetric polynomial of degree at most iV(iV — l)/2. 
It must be divisible by the Vandermonde determinant Vjy(X), and it remains to 
verify the constant prefactor. Following the highest in lexicographic order term 
Xi~ l X2~ 2 ■ ■ ■ x° N we see that upon the action of (Pi + ■ ■ • + T>n) it collects the 
coefficient 

JV-l 

^2{jij-i)+sj), 

which sums to (15. lOj) . Thus, (15. lip is proved. 
Step 2. Let us now prove that 

p Nit; A, i/) = 1, A G GTjv, t>0, (5.17) 

with P N as in (15.121) . In the space Xjv, (I5.12p reads 

P N (t;X,Y) = e-^^^detlP^x^)]" X, Y e X N . (5.18) 

Since the action of T> in the space of polynomials lSL[x] is consistent with filtration 
by degree, see (I5.16p . the action of the corresponding semigroup (-Pi(^))t>o in M[x] 
is well-defined, and (I5.16P implies 

x, y)y m = e {m{m ~ l)+sm)t x m + lower degree terms. 

y&1 
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We obtain 



= vTx) z\2 sgna X, p i(*;^(i)'2/i) ■ ■ • p i(t'i x <T<.N),VN)v?~ 1 v£~ 2 - ■ -un- 



e~ dNt \ - ,^.v 



£ sgn.e^^O'O'-i)^)^^ . . . Xct(7V _ 1} = 1; (5 . 19) 



where Sn denotes the group of permutations of {1, ... , N}. Note that the first equal- 
ity (change of the summation domain) holds because the expression for P/v(t; X, Y) 
is symmetric in (Uj), and it vanishes if yi = yj for i ^ j. 

Step 3. Consider N independent copies of the bilateral birth and death process of 
Subsection 15. 1[ and denote by 7i n (t;X, Y), 1,7 6 3C N , the probability that these 
processes started at x\, . . . , xjy end up at y%, . . . , after time t having made a total 
of n jumps all together, and their trajectories had no common points at any time 
moment between and t. We want to show that 

P% ] (t;X,Y) = e-^^7r„(t;X,n (5.20) 

where Pff is defined as in Section H] using P^' as the matrix of transition rates. 
Indeed, computing 7r n 's boils down to recurrence relations 

n (t;X,Y) = e^^ Y h x=Y , 
; X, Y) = f e ^l"W) Y, V inl( X > Z >n-i{t - r; Z, Y)dr, n>l, 

where 'D^J = T>\ + • ■ ■ + T>n is the matrix of transition rates for the N independent 
birth and death processes. 

For n = 0, (15.201) follows from ( 15. 9p . Assuming ( I5.20p holds for n — 1, we rewrite 
the recurrence relation for 7r ra 's as 

7r n (t;X,Y)= fe T(2,W( ^W 
Jo 

Comparing with the recurrence relation for P^ n \ cf. Section HI yields (I5.20p . 



MARKOV PROCESSES 



31 



Step 4- Following Section H] and using ( 15.20(1 . we see that the minimal solution for 
the backward equation with T>^ N ' as the matrix of transition rates, has the form 

Mt; X,Y) = jr P^ ] (t; X, Y) = e^^S £ n n (t; X, Y). 

n=0 ^ ' n=0 

The last sum is clearly equal to the probability that N independent copies of the 
bilateral birth and death process of Subsection 15.11 started at xi,...,xn end up at 
yi, . . . , after time t without intermediate coincidences and without any restriction 
on the number of jumps. Note that we are using the fact that the birth and death 
process does not make infinitely many jumps in finite time (minimal solution is 
stochastic), cf. Theorem 15. II 

Such a probability of having nonintersecting paths is given by a celebrated formula 
of Karlin-McGregor (KM59j : 

oo 

7r n (t; X, Y) = det [P^t; x h Vj j\ £. =1 , X, Y, G X N . 

n=0 

Hence, the minimal solution P/v(t; X, Y) coincides with the right-hand side of (15.181) . 
and by Step 2 it is stochastic. We have thus shown that the matrix of transition 
rates on GTjv is regular, and the semigroup has the form ( 15.121) (or (15 . 18(1 ) . 

Step 5. To conclude the proof of Theorem 15.41 it remains to show that the Markov 
semigroup (-P/v(t))i>o is Feller. This is equivalent to proving that 

lim P N (t; A, u) = 0, t > 0, v G GUV (5.22) 

A— >oo 

But this immediately follows from ( 15. 12(1 because we already know that ( 15.221) holds 
for N = 1 (Theorem 15. ip . and DimTv(A) is always at least 1. 

□ 

6. COMMUTATIVITY 

The goal of this section is to address the question of compatibility of the semi- 
groups of Section [5] and links of Section [31 cf . (12.31) . 

6.1. Parameterization. As we shall see, in order for the commutativity relations 
(12. 3p to be satisfied, the parameters (u, u', v, v') used to define semigroups (Pjv(£))t>o 
need to depend on N . For that reason, introduce two new pairs of parameters (z, z') 
and (w,w r ) that satisfy the same conditions as (u, u', v,v') before: 

(z + k){z' + k) > 0, (w + k)(w' + k) > VA; G Z; z + z' + w + w' > -1. 

(6.1) 

Furthermore, for N > 1 define 

u N = z + N-l, u' N = z' + N - 1, v = w, v' = w', (6.2) 
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and let (-Pzv(£))t>o be the Feller semigroup of the previous section with parameters 
(u,u',v,v') = (u N ,u' N ,v N ,v' N ). 
We are aiming to prove the following statement. 

Theorem 6.1. With links {A% +1 } N>1 as in Subsection\3J\ and semigroups (-P/v(£))t>o 
as above, the compatibility relations ( 12. 3 p hold. 

6.2. Infinitesimal commutativity. We first prove a version of (12.31) that involves 
matrices of transition rates. 

Proposition 6.2. For any N > 1, u,u',v,v' G C, and A G GlV+i,^ G GTn, we 
have 

& N+r> (A, k)A^ +1 (k, u) — Af(A,^ W (^) (6.3) 
or, in matrix notation, 

V (n+i) a n+i = A n+i v (n) ; where V (N) is the opemtor defined 

by (15.91) . and in p( 7V+1 ) we replace N by N+l and the parameters (u, u') by (u, u') = 
(u + l,u' + 1). 

Proof. We start with the following simple lemma. 

Lemma 6.3. Let [A (A, u)\ AgGTjv+i „ gCTjv be a nnatrix with rows parameterized by 
GT^v+i and columns parameterized by GUV, and such that each row of A has finitely 
many nonzero entries. If for any symmetric polynomial F in N variables and any 
A G GTat + i we have 

A{X,u)F{v L + N-l,u 2 + N-2,...,u N ) = Q, (6.4) 

then A(X, v) = 0. 

Proof. Assume A(\, v) ^ for some A and v. Let z^ 1 ) , . . . , z/^ G GTat be all 
signatures different from z> and such that A(\, u^>) ^ 0. 
Set x = {j>x + N - 1, . . . , i) N ) G Z N and 

y ^=(^ + N-l,...,u^)eZ N , j = l,...,L 

Observe that the orbits of the vectors x,y^\ . . . ,y( l > under the group of permuta- 
tions of the coordinates do not intersect. It follows that there exists a polynomial 
/ in N variables, which takes value 1 on the orbit of x and vanishes on the orbits 
of the vectors . . . , y®. Then for the symmetrized polynomial F(zi, . . . , zn) = 
J2aes N fi z c{i)i ■ ■ i z a(N)) the left-hand side of (16 .4p is equal to N\A(\, v) ^ 0. Con- 
tradiction. □ 
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Let us now introduce symmetric polynomials on which we will evaluate (in the 
sense of Lemma both sides of (16.31) . For a partition (=signature with nonnega- 
tive coordinates) \i G GT N and c G C set 



(-^")yLt \ \l <~~ ^ AT {Xj X -j 



l<i<j<N\^i ^3 

1 d € 



1 detr(x, + c)^ +Ar+l -^l 7V + l 1 



(n + i) M n^i^xTv+ii^ - x i) 

where we assume Hn+i — and use the notation (a G C, G Z> ) 
a^ fc = a(a — 1) ■ ■ ■ (a — fc + 1), (a) k = a(a + 1) • • • (a + k — 1), a^ = (a) = 1, 

A 7 " 

(a) M = JJ(a-j + l)^. 

i=i 

Clearly, and (j MiC are symmetric polynomials in N and N + 1 variables, respec- 
tively. Moreover, for any fixed c G C, the polynomials {F^^} with /i ranging over 
all nonnegative signatures in GTpj form a linear basis in the space of all symmetric 
polynomials in N variables. Indeed, this follows from the fact that the highest degree 
homogeneous component of F^ c coincides with the Schur polynomial s^(xi, . . . ,Xjy), 
and those are well known to form a basis, see e.g. |Macd95j . 

Hence, to prove Proposition [672] it suffices to verify that the two sides of ( 16. 3ft give 
the same results when applied to F^ c for a fixed c and // varying over nonnegative 
signatures of length N. 

Lemma 6.4. For any A G GTjv+i, any nonnegative signature \i G GT^, and c G C, 
we have 

A^ +1 (A, v)F fttC (u 1 +N-l,...,u N ) = G^ c {\ x + N, A 2 + N - 1, . . . , X N+1 ). 

v£&T N 

Proof. The argument is similar to that for relation (10.30) in [0097J. Denote 

(xi, . . . , x N+1 ) = (Ai + N, . . . , Ajv+i), (yi, ■ ■ ■ , Vn) = (vi + iV - 1, . . . , v N ). 

Then z/ -< A means < ?/j < Xj for all i = 1, . . . , N. Taking into account the 
definition of A^ +1 , one sees that the relation in question is equivalent to the following 
one 

det\( Xi + c)^ +N+1 ^] N+ ^ 

= £ det[( K + c)«««-»]^ r (6.5) 

1 yi ,...,y N & 

Xi^.i<yi<Xi for all i 
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The last column in the (N + 1) x (N + 1) matrix in the left-hand side of ( 16. 5ft 
consists of l's. Subtracting from the zth row the (i + l)st one for each i — 1, . . . , N, 
we see that the left-hand side is equal to the N x N determinant 



det 



{xi + c + _ + c+ ^Um+N+l-j) 



N 



»J=1 



On the other hand, the summation in the right-hand side of ( 16. 5 j) can be performed 
in each row separately using the relation 

( 6 + c )4(m+l) _ ( fl + c )4(m+l) 



m + 1 

y=a 

Collecting constant prefactors completes the proof of Lemma [6 

(iy + iyv! 

= MT] { " J • V • ' ■ /S!( V - /S! = 1. (6.6) 

□ 

To conclude the proof of Proposition 16.21 we want to prove that, for a suitable 
fixed constant c G C, T>^F^ C decompose on {F uc \ in exactly the same way as 
V^ N+1 ^G^ C decompose on {G u<c }. 

It is actually convenient to take c = v , where v is one of the four parameters 
(u,u',v,v f ). With this specialization we prove 

Lemma 6.5. For any X G GT/v and any nonnegative signature n G GTjv, with the 
notation rrij = fij + N — j , j = 1, . . . , N , we have 

V ( ^ N \X,u)F ll:V (u 1 + N-l,... } iy N ) 

u€GT N 

N N \ 

^ rrij (rrij - 1) + s^rrij - d N j F^(Ai + iV - l,...,Ajv) 

vi=i j=i / 

JV 

+ ^ ({rrij - l)(v' - v + rrij - 1) + s(rrij - v - 1) + uv! - vv'j l H -i> H+1 

xF At _ e ., 1> (A 1 + JV-l,...,A JV ), (6.7) 

where is as in (j5.10j) . = (0, . . . , 0, 1, 0, . . . , 0) mi/j 1 at £/ie jt/i place, and we 
assume (jln+i = 0. 
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Proof. We first compute, cf. (15.1 5ft . 



35 



V(x + v) lm = (x + v)(x + i/)AV(i + v) im + (sx + uu' - vv')A(y + v) im 

= m(m -l)(x + v')(x + v) i(m ~ 1) + m{sx + uu' - vv') (x + v)^ m - 1} . (6.8) 

This is the place where the choice of c = v matters; for different values of c the 
expression for T>(x + c)^ m would have been more complicated. 
Substituting 

x + v' = (x + v — m + 1) + (v' — v + m — 1), 
sx + uu' — vv' = s(x + v — m + 1) + (s(m — u — 1) + uu — vv'), 

we obtain 

V(x + v) lm = (m(m - 1) + sra) (x + v) 4 ™ 

+ ((m - l)(w' - d + m - 1) + s(m - v - 1) + uu' - vv')m(x + v)^ m ~ 1 \ (6.9) 

The statement now follows from (15.91) and the definition of F^ c . 

□ 

Let us complete the proof of Proposition 16.21 

Apply both sides of ( 16. 3D to F^ v in the sense of Lemma 16.31 Using Lemma 16.41 
we see that the left-hand side of (16.31) turns into 

& N+l \x, ^G^im + n,..., K n+1 ), 

and repeating the arguments of Lemma 16.51 we see that this is equal to 

' N N \ 

Y rhjirhj - 1) + s ^ rhj - d N+1 J G^Xi + N - 1, . . . , X N ) 

N 

+ ^2 W^i - l){v' -v + rhj - 1) + s(rhj - v - 1) + uu' - vv'j l w _i> w+1 

3=1 

xG^ Sj , v (Xi+N-l,...,X N ), (6.10) 

where rhj = Uj + N +1 — 1 = rrij + 1, and tildes over the other constants mean that 
in their definitions we replace (u, u') by (u, u') = (u + 1, u' + 1). 

On the other hand, by Lemmas 16.41 and 16.51 the right-hand side of (16.31) equals 
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N N \ 

22 m j( m j _ 1) + s '^2 m j ~ I G^ v {\i + N - 1, . . . , Aat) 
0=1 3=1 J 

N 

+ (jjn>j — l)(v' — v + rrij — 1) + s(rrij — v — 1) + mm' — ff ' 

3=1 

x l^_i>^ +1 <V^( A i + iV - 1, . . . , Aiv). (6.11) 

It is a straightforward computation to see that all the coefficients in ( 16.101) and 
(16. lip coincide. The proof of Proposition 16.21 is complete. 

□ 

6.3. From matrices of transition rates to semigroups. In order to complete 
the proof of Theorem 16.11 we need the following lemma. 

Lemma 6.6. If f is a finitely supported function on G1V then A^ +1 f is in the 
domain of the generator An +1 of the semigroup (P/v + i(t))i>o (see Section^ for the 
definition of the generator and its domain). 

Let us postpone the proof of Lemma 16.61 until the end of this subsection and 
proceed with the proof of Theorem 16.11 

In order to prove (I2.3p it suffices to prove that the two sides are equal when applied 
to a function / on GT N with finite support (as such are dense in C (GT N )): 

P N+1 (t)A^ +1 f = A% +1 P N (t)f, t>0, iV = 1,2,.... (6.12) 

Let us denote the left and right-hand sides of (16.121) by Fi e f t (t) and F rig ht(t). We will 
show that they solve the same Cauchy problem in the Banach space Co(GTtv+i). 
Then (16.121) will follow from an abstract uniqueness theorem for solutions of the 
Cauchy problem for vector functions with values in a Banach space, 

j t F(t) = AF(t), t > 0, F(0) = fixed vector, 

which holds under the assumptions that (1) A is a closed dissipative operator, (2) 
F(t) is continuous for t > and strongly differentiable for t > 0, and (3) F(t) G D(A) 
for t > 0; see e.g. |Kat80l IX. 1.3]. 

In our situation, A = A^+i and the fixed vector is A^ +1 /. Obviously, both Fi e f t (t) 
and F rig ht(t) are continuous for t > and they have the same initial value A^ +1 / 
at t = 0. 

Let us check the differential equation for Fi e f t {t). By Lemma IBTBI we have A^ +1 / £ 
D(A N+ i). Hence, F [e f t (t) £ D(A N+ i) (semigroups preserve the domains of the 
generators) and it satisfies 

j t Fieft(t) = A N+1 F left (t), t > 0. 
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Let us turn to F rig ht- By Corollary 14. 7\ f belongs to D(An). It follows that the 
function t h-> Pjv(t)/ is strongly differentiable and 

j f PN(t)f = A N P N (t)f = vWp N (t), t > 0. 
Hence, F right {t) is also strongly differentiable and for t > 

jF nght {t) = ^ N +l jPN{t)f = A% +1 vWp N (t)f. 

By definition, the last expression should be understood as A^ +1 (V^ N \P N (t)f)). 
However, since all rows of the matrices A^ +1 and T>^ N ' have finitely many nonzero 
entries, we may write 

K N N +l {V^\P N {t)f)) = (A»+ l vW)P N (t)f. 

By virtue of Proposition 16.21 this equals 

V^A^P N (t)f = & N +^F right (t), 

so that 

-j^F r i g ht(t) = U( N+1 ^ F rig ht(t). 
Next, as ^F r i g ht(t) is in Co(GT/v+i), so is 

V (N+1) 

F r ight(t). By Proposition |461 we 
may replace T>( N+1 ) by A^+i, which gives the desired differential equation 

-j^F right (t) = A N+ iF right (t), t > 0, 

and we conclude that F ie f t = F right . 

Thus, we have proved Theorem 16.11 modulo Lemma [6.61 

Proof of Lemma 1 6. b\ Let / be a finitely supported function on GTtv, g = A^ +1 f. 
Proposition 13.31 says that g G C (GT N+ i), and by Proposition 14. 61 it suffices to check 
that V( N+1 *>g G C (GT N+1 ). We have 

(V^g)(X) = Yl ^ (W+1) (A, A + e) (g(X + e) - g(X)), 

e:X+e£&T N+1 

where A G GTjy+i, £ ranges over {±ej}j = i ) ... ) jv+i, with (e,) being the standard basis 
in IR 7V+1 , and T)( N+1 \\, A + e) are off-diagonal entries of the matrix T>( N+1 \ 

Without loss of generality we may assume that / is the delta-function at some 
v G CIV. We obtain 

t\\- m ni<i<j<ivO; -i-v,+ j) 

g{* — -Ff 7T : 7 — — ~v ■ K^X, 

lll<i<j<N+l( X i - 1 - X j+J) 
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and 

where = A 3 - + iV + 1 — j, j = 1, . . . , N + 1, and we assume Dimjv+i(A + 6^) = 
in case A + 5^ ^ GTtv+i (this is supported by the explicit formula for Dirn/v+i( • )). 
Observe that for 



HA) 



Nl Hi<i<j<N( u i j) _ consti 



rii<i<j<JV+i( A i - i - X j + i) Dim w+1 (A) 
(we removed the factor l u ^\ from g(X) above), we have 



const 5 



Dim 7v+i(A) 

V(k,k + ei)(Dim N+l (\ + e iei ) -Dim w+1 (A)) = Drr^lT) ' (6 ' 14) 



x 

i=\ e»=±l 



where we used (I5.14p . 

Next, observe that the function 



(^ +1 ^)(A)W = D . c ° nst3 rA . W 

DmiAr+i(A) 

belongs to Co(GTjv+i). Indeed, if A goes to infinity inside the subset {A : v -< A} then 
Aj— Xj — > +oo for at least one couple i < j of indices, which entails DimAr + i A — > +oo. 

The discrepancy between (V^ N+1 ^ ) 'g)(X)l u ^x and (T'( 7V+1 ^)(A) (or rather between 
the summations in ( 16. 131) and (16. 14ft ) comes from values of i and £j such that either 
v -< A but z/ ^ (A + e^ei), or z> -< A + £jej but v -ft X. In both cases, for that value of 
i, the quantities Aj, h, and /j + £,) must remain bounded as z/ is fixed. 

Note that A — > oo inside the subset 

{A G GTjv+i :z/-<Aorz/-<A + £jej for some z}, 

then either Ai — >■ +oo or Xn+i — > — oo, or both, while all other Xj remain bounded 
from both sides. But then a direct inspection of the summands in ( 16 . 13f) and (16. 14ft 
that contribute to the discrepancy shows that they converge to zero as A — > oo. 
Hence, (V^ N+1 ^g)(X) G C (GT N+1 ). * ' □ 

7. Invariant measures 

In three previous sections we defined a chain of countable sets {E N = GT N } N >i, 
constructed links A^ +1 between then, and identified the boundary = f2. Fur- 
thermore, for any quadruple of complex parameters (z, z' ,w,w') satisfying (16. ip we 
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constructed Feller semigroups (-P/v(£))t>o on GTjy and showed that they are com- 
patible with the links; by Proposition 12.41 this yields a Feller semigroup (P(t)) t > on 

a 

The goal of this section is to exhibit an invariant measure for (P(t)) t >o- 

7.1. ^-measures. Let z,z',w,w' be complex parameters satisfying (16 . 1 p . As was 
pointed out in Subsection I5.1[ this is equivalent to saying that each pair (z, z') and 
(w, w') belongs to one (or both) of the sets 

{(C,Oe(C\Z) 2 |C' = C} and 

{(C, CO e (K \ 1 m < C, C' < m + 1 for some meZ}, 

and also z + z' + w + w' > —1. 
For XeGT N set 

M ZtZ i >w>w i\ N (\) = (const^) -1 • M' Z)Z , yWjW ,\ N (\) 

where 

n ( 1 
m 'z,z',w,w>\n{X) = II [ T{z _x i + i)T(z'-X i + i) 

r(w + n + 1 + \i — o!v + n + 1 + a, - i) ) • ( Dim ^ A )) 2 ' (") 

and 

const w = ^ M L',™xi^( A ) 

AeGTjv 

is the normalizing constant depending on z, z', w, w', N. 

Theorem 7.1 QOls03j). Under our assumptions on the parameters, for any N > 1, 
M z>z t w w iw is a probability measure, we call it the Nth zw-measure. Moreover, these 
measures are consistent with the links, 

^-z,z' .w,w'\N = M z ^ z i )W:W i\ N+ i A^ +1 , N > 1, 

with A^ +1 as in Subsection \3. 1[ 

Theorem 17. II implies that the system {M ZjZ i )WtW n^)^>i defines a probability mea- 
sure M ZiZ ' )W}W i on the boundary Q that we call the spectral zw-measure, cf. Theorem 
13.1^ and a character of the infinite-dimensional unitary group U(oo), cf.[Ols03]. For 
z' — z and w' = w one can find a geometric construction of the corresponding 
representations of U(oo) in [Ols03] . There is also a fairly simple "coordinate-free" 
description of general zw-measures that we now give, cf. [B O05c| . 

Let T be the unit circle in C and be the product of iV copies of T (the 
A^-dimensional torus). For any A G GTjy, the character x X of the corresponding 
irreducible representation tt\ of U(N) can be viewed as a symmetric function on T , 
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where coordinates are interpreted as eigenvalues of unitary matrices. Explicitly, the 
character is given by the (rational) Schur function 



X (lii,...,Utf) = S X {Ui,...,U N ) - 



Xi+N-ji 

det[«f^] : 



l<i,j<N 

Consider the Hilbert space of symmetric functions on T , square integrable 
with respect to the measure 

1 n 

Wl n \^-^\ 2 u du 



i ■ 



l<i<j<N i=l 



which is the push-forward of the normalized Haar measure on U(N) under the 
correspondence U h-> (u±, . . . ,un). Here dui is the normalized invariant measure on 
the ith copy of T. 

Given two complex numbers z, w, we define a symmetric function on by 



N 



fz,w\N{ u ) = Yli 1 + Ui) Z (l + Ui 



i=l 



If $l(z+w) > — ~ then f z>w \N belongs to the space H N . Let (V, w') be another couple 
of complex numbers with + w') > — \. We set 



(fz,w\N, Xa)(XA) fo,V\NJ 



M Z>Z / >W7W >\ N (\) — -r- 't , A G GTjy, 



where ( • , • ) is the inner product in H N . It turns out that this definition leads us 
to the explicit formula given above. 

The spectral zw-measures were the subject of an extensive investigation in [BO05aj 
the upshot of which is the statement that with u e Q distributed according to 
M ZtZ i tW ^ w i, its coordinates 

(where possible zero values of af and /if should be removed) form a determinantal 
point process on R\ {±|} with an explicit correlation kernel. See |BO05aj . |BO0 5b] 
for details. 

7.2. Invariance. The main statement of this section is 

Theorem 7.2. For any quadruple (z, z',w,w') of parameters satisfying (16. ip . the 
spectral zw -measure M ZtZ t >WyW > is the unique invariant probability measure with respect 
to the semigroup (P(t)) t >o- 
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Proof. Let us prove the invariance first. By Subsection 12.81 it suffices to verify that 
for each N > 1, the iVth level zw- measure is invariant with respect to (P/v(^))i>o- 
We will check this fact on the level of matrices of transition rates: 

E M Z;Z , )W>w , lN {X)V( N \\,v) = 0, N>1, ueGT n . 

AgGTjv 

Since it is easy to check that P^' is reversible with respect to M z>z i jW>v} i\n(X) } 
M Stg ,^\ N {\) D (A °(A, v) = V (N \\, v) M z>z , tW , w , lN (iy), X,ue GT N , 

an argument in Section 3 of [Kel83] shows that the invariance on the level of tran- 
sition rates implies the invariance with respect to the corresponding semigroup. 

As in the proof of Theorem I5.4[ it is convenient to employ the bijection A <H- 
(Aj + N — j)i<j<N between GT N and Xn, see Subsection 15.21 for the notation and 
also recall that we are using parameterization (16.21) . Under the bijection of GTjv 
and Xat, the desired identity takes the form (removing irrelevant prefactors) 



e ( n ) (w^ + • • • 

h V(x N , y N )l {Xi=m ^ N} ) - d N l x =Y ) = 0, (7.2) 



X£X N \i=l 



where 

v; T(z + N - x)T(z' + iV - x)T(w + 1 + x)T(w' + 1 + x) ' 

Let po = 1 , pi , P2 1 ■ ■ ■ , deg pj — j — 1 , be monic orthogonal polynomials on Z 
corresponding to the weight function W(x). As 

W(x) = 0(\x\- z ~ z '- w - w '- 2N ), x^oo, 

the assumption z + z' + w + w' > —1 implies that W(x) has at least 2N — 3 finite 
moments, and polynomials pj with j = 0, 1 . . . , N — 1 are well defined. 

Polynomials {pj} can be written explicitly in terms of the hypergeometric function 
3_F 2 evaluated at 1. They were discovered by R. Askey [Ask87j . and independently 
by P. Lesky |Les97j . |Les98j ; see also the recent book |KLS10t §5.3, Theorem 5.2, 
Case IIIc]. We call them the Askey-Lesky polynomials . 

The Askey-Lesky polynomials are eigenfunctions of the operator T> on Z, see 
[BU05al §7]: 

E v)Pj(y) = ^jPj( x ) Vx g z, j = o, i, 2, . . . , 

yei, 

where 

Ij = j(U - 1) - (un + u' n + v n + v' n )). 
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Multiplying both sides by W(x) and using the fact that W(x)V(x, y) is symmetric 
with respect to transposition x 4-> y we obtain 

Y,Pj{xW{x)V{x : y)= ljPj {y) VyeZ, j =0,1,2,.... (7.3) 

Let us rewrite the Vandermonde determinant in the left-hand side of (17. 2p as 

V N (x) = ±dek[p i - 1 (x j )] itj=1 . 

Applying operators T>i, . . . ,T>n to individual columns in this determinant multi- 
plied by W(xi) ■ ■ ■ W(xn) according to (17.31) . and recalling the definition of d N , we 
obtain (17211 . 

Let us now prove uniqueness. As explained in Subsection I2.8[ it suffices to 
show that the iVth level zu>-measure is the unique invariant probability measure 
for (-P/v(£))t>o f° r an Y N > 1. But uniqueness of invariant measures holds in general 
for irreducible Markov chains on countable sets, see e.g. Theorem 1.6 in |And91j . 

□ 

8. Stochastic dynamics on paths. General formalism 

8.1. Overview. Let us return to the general setting of Section [2] and assume that 
all -EVs are discrete. For N — 1,2, . . . set 

N-l 

E (N) = [(xt, . . . , x n ) e Ex x • • ■ x E N | 11 At + \x k+1 , x k ) ± o}. (8.1) 

k=l 

There are natural projections n^ +1 : E^ N+1 ^ -»■ EW consisting in forgetting the 
last coordinate; let E^°°' = hmE^ N ', where the projective limit is taken with respect 
to these projections. Obviously, E^°°> is a closed subset of the infinite product space 
HivLi En- Thus, elements of E^°°' are some infinite sequences. Let LT^ : E^°°' — > 
EW be the map that extracts the first N members of such a sequence. 

Definition 8.1. We say that a probability measure on E^ N ' is central if there 
exists a probability measure [In on E N such that 

/i (7V) (xi, . . . , X N ) = /iAr(x A r)A^_ 1 (x A r, X N -l) ' ' ' ^l(x 2 , X X ) (8.2) 

for any (x\, . . . ,xn) G E^ n \ Relation (I8.2p establishes a bijection between proba- 
bility measures on En and central probability measures on 

We say that G M. P (E^ 00 ' > ) is central if all its pushforwards under projections 
are central. Relation ( 18. 2 p also establishes a bijection between central measures 
on E^ and elements of \^mAi p (E N ) of Subsection 12.51 

Finally, we say that a Markov semigroup (P( N \t)) t >o on E^ is central if the 
associate linear operators in Jvl^E^) map central measures to central measures. 
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Clearly, a central Markov semigroup (P^ (t))t>o defines a Markov semigroup on 
E N — in order to obtain jx N P N {t) for fx N G M. p (E N ) one needs to define via 
dS2|l, evaluate /iWpW(f), and read off a measure on £jv using Definition 18.11 

Proposition 8.2. Le£ (P (A °(t)) <>o? N > 1, be a sequence of central Markov semi- 
groups on E( N ^ 's that are compatible with the system of projections: 

p( N+1 \t) o n^ +1 = n^ +1 o p( N \t), t > o, n > 1. 

Then the corresponding Markov semigroups (-P/v(t))t>o on En, N > 1, are compat- 
ible with projections A^ +1 as in (I2.3p . 

Proof. Follows from the fact that if iv- N+l > and hn+i are related as in Definition 18. II 
then yu( iV+1 )ll^ +1 and /iAr + iA^ +1 are also related in the same way. □ 

The goal of this section and the next one is to construct central Markov semigroups 
(P {N) (t)) t >o that would yield, as in Proposition 18.21 semigroups (-P/v(t))t>o on En = 
GTn that we dealt with in the previous sections. One reason for such a construction 
is the fact that for the Gelfand-Tsetlin graph, the isomorphism between central 
measures on GT' 00 *' and probability measures on the boundary Q, cf. Definition 12.21 
is somewhat explicit, see Section below. Thus, (P( N \t)) t > can be thought of as 
providing a more "hands-on" description of the corresponding semigroup (P(t)) t > 
on Q. 

8.2. Construction of bivariate Markov chains. Let E and E* be countable 
sets, and let Q and Q* be matrices of transition rates on these sets. Let A = 
[A(x*, x)] x *ee*,x<ee be an additional stochastic matrix which we view as a stochastic 
link between E* and E. 

We will assume that for each of the three matrices Q, Q*, and A, each row contains 
only finitely many nonzero entries. In addition, we assume the relation 

Y,A(x*,x)Q(x,y) = ^ Q*(x*, y*)A(y*, y), x* eE*, ye E, (8.3) 

xeE y*eE* 

or AQ = Q*A in matrix notation. 

Observe that in case A(x*,y) = 0, the diagonal entries Q(x,x) and Q(x*,x*) give 
no contribution to (18. 3p . and the commutativity relation can be rewritten as 

A(x*, x)Q(x, y)= Q*( x *i y*) A (y*, V), x* G E*, yeE. (8.4) 

x£E,x^£y y*£E* ,y*^x* 

We will denote the above expression by A(x*,y); it is only defined if A(x*,y) = 0. 
In what follows we also use the notation 

qx = —Q(x,x), x G E; g*, = — Q*(x*, x*), x* G E*. 

Consider the bivariate state space 

£( 2 ) = {(x*,x) <E E* x E \ A(x*,x) ^ 0}. 
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We want to construct a Markov chain on that would satisfy two conditions: 

• The projection of this Markov chain to E gives the Markov chain defined by Q; 

• It preserves the class of measures on E^ satisfying Prob(x|x*) = A (a;*, x)\ 

• In this class of measures, the projection of this Markov chain to E* gives the 
Markov chain defined by Q*. 

To this end, define a matrix Qw of transition rates on E^ with off-diagonal 
entries given by 



Q®((x* 7 x),(y*,y)) 



'Q(x,y), 
Q*(x%y 



Q{x,y) 



^ My*,*) 

A(x*, x) ' 
Q*(x*,y*)A(y*,y) 

A(x*,y) 



,0, 



x' = y, 

x = y, 

A(x*,y) = 0,A(x*,y)^0, 
otherwise. 



Note that A(x*,y) = implies x* ^ y* and x ^ y (provided that (x*,x), (y*,y) 
are m E^) so all the cases in the above definition are mutually exclusive. 
The diagonal entries Q^((x*,x), (x*,x)) with (x*,x) G E^ are defined by 



-Q^{(x*,x),(x*,x))=q£\ x) := ^ ((**, x), (y\ y)) . 

Clearly, any row of also has only finitely many nonzero entries. One immediately 
verifies that for any (x*, x) G E^ and y G E with x ^ y, 



^ QW((x*,x),(y*,y))=Q(x,y). 



i.5) 



Indeed, one needs to consider two cases A(x*,y) = and ^ 0, and in both cases the 
statement follows from the definitions. As the row sums of and Q are all zero, 
we obtain (I8.5P for x = y as well. 

For any x G E, let us also introduce a matrix of transition rates Q x on the fiber 
E x = {x* G E* | A(x*,x) ^ 0} via 



Q x (x*,y*) = Q {2 \(x*,x), (y*,x)) = Q*(x*,y*) 



*n A(y*,x) 



A(x*,x) ' 



x t y > 



and 



y*eE x ,y*^x* 

The following statement is similar to Lemma 2.1 of |B F08+ | proved in the dis- 
crete time setting. As we will see, the proof of the continuous time statement is 
significantly more difficult. 
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Proposition 8.3. Assume that the matrices of transition rates Q, Q* , and Q x for 
any x G E are regular. Then is also regular, and denoting by P(t), P*(t), and 
p( 2 \t) the transition matrices corresponding to Q,Q*, and Q^ 2 \ we have 

P [2) (*l « x), {y*, y)) =P(t- x, y), (8.6) 

£ A (^*' x ) p{2) (*; *)> = p *(*; **> ( 8 - 7 ) 

z:(x*,z)G£( 2 ) 

where in the first relation (x*,x) G E 1 ' 2 ', y £ E are arbitrary, while in the second 
relation x* G E*, (y*,y) G E^ are arbitrary. 

Proof. The regularity of and collapsibility relation ( 18. 6ft follow from Proposition 
S3 with (EC2) specializing to ( 1831) . 

Proving (18 .7p is more difficult, and we will follow the following path. First, we will 
show that both sides of (18. 7p satisfy the same differential equation (essentially the 
Kolmogorov backward equation for P* (£; x* , y*)) with a certain initial condition. 
Then we will see that the right-hand side of (18 ,7p represents the minimal of all 
nonnegative solutions of this equation. Since for a fixed x*, both sides of (18. 7p 
represent probability measures on E^ 2 \ the equality will immediately follow. 

For the first step, let us show that the left-hand side f t (x*,y*,y) of (18 ,7p satisfies 

j t ft{x\y\ y) = Yl Q*( x *> **)/*(**, V*. v) (8-8) 
z*eE* 



with the initial condition 



l\m Q f t (x*,y*,y) = l x * =r A{y*,y). (8.9) 



The initial condition satisfied by P^ 2 '(t) implies ( I8.9p . so let us prove 
Using the Kolmogorov backward equation for p( 2 \t), we obtain 



d 
dl 



ft(x*,y*,y) = £ A(x*,x)(-q^ )X) P^(t;(x*,x),(y*,y)) 

x:(x*,x)£E( 2 ) \ 

+ Yl Q {2 \tx\x)i{z\z))pV\t-{z\z),{y\y))\ (8.10) 

(z* ,z)Mx* ,x) / 



(z*,z)=£(x*,x) 

For the first term in the right-hand side, we use 

A(w*, x 



A(x*,x) 

w*:w*^x* 
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which follows directly from the definition of Q^. Thus, we can rewrite the first 
term in the right-hand side of (18.101) as 

-J2<lxMx*,x)pW{t; (x*,x),(y*,y)) 

x£E 

- E Hw*,x)Q*(x\w*)P (2) {t;(x*,x),(y*,y)). (8.11) 

x: (x*,x)£E( 2 ) w*:w*jtx* 

For the second term of the right-hand side of (18.101) . according to the definition 
of Q( 2 \ let us split the sum over (z*, z) into three disjoint parts: (1) x* = z*, x ^ z; 
(2) x* ^ z*,x = z; (3) A(x*, z) = (hence, x* ^ z*, x ^ z). 

Part (1) gives 

(1) = A ( x » E Q(x,z)pW{t;(x*,z),(y*,y)). 

i:(i',i)eE( 2 ) 2:2/i,(i',z)6E( 2 ) 

Interchanging the summations over x and z, we can employ the commutativity 
relation (18.31) . This gives 

(1) = E E A « z)Q*(x*,v*)P {2) (* ; (x*, z), (y*, y)) 

2:(i*, 2 )eB( 2 ) v*£E* 

+ vMx*, x)P {2) (t; (x*, x), {y\ y)) . (8.12) 

x£E 

Observe that the last term cancels out with the first term in (18.111) . while the sum 
of the first term of H8. 121) and the second term of (18. lip , with identification z — x, 
v* = w* of the summation variables, yields (only terms with v* = x* survive) 

- «£. H**, x)P {2) (t; (x*,x), (y*, y)) . (8.13) 

xGE 

Further, part (2) of the second term of ( 18. lOf) reads 

(2) = E E Mz*,x)Q*(x*,z*)pM(t;(z*,x) } (y\y)). (8.14) 

i:(i*,i)eE( 2 ) z*:z*^x* 

Finally, part (3) gives 

(3)= E E h^m*,*) q ' {x '£T z) 

x:(x* ,x)&E( 2 ) {z* ,z)£EW :A{x* ,z)=Q 1 ' ' 

xPW(t; (z*,z),(y*,y)) 
= E E Hz*,z)Q*(x*,z*)P (2) {t;(z*,z),(y*,y)), (8.15) 

z:(x* ,z)<£E( 2 ) z*:z*^x* 
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where we used the definition of A, see ( 18.41) . to perform the summation over x 7^ z. 
One readily sees that adding f)8.13j) . f)8. 14j) . (18. 15j) yields the right-hand side of ( 18. 8p . 

Assume now that we have a nonnegative solution f t (x*,y*,y) of ( 18. 8 p satisfying 
the initial condition ( 18. 9p . Multiplying both sides of ( 18. 8 p by exp(q* t t) we obtain 

(exp(q:j)f t (x\y\y)y = exp(^i) £ Q\x\z*)f t {z\y\y). 
Integrating both sides over t and using (18. 9p gives 



ft(x*,y*,y) = l x * =y *A(y*,y)exp(-q**t) 

+ f exp(-q* x *s) V Q*(x*,z*)ft- S (z*,y*,y)ds. (8.16) 
Jo z « +x * 

Set F} \x*, y*) = l x * =y * exp(— g*,t), and for n = 1, 2, . . . define 

f/ b V,v*) = ^ ( V,v*) + teM-q» E < 5*( x *' z *) F t- s ~ 1) ( 2; *^*) rfs - 

Clearly, (18. 16|) implies f t (x*,y*,y) > Ft (x*, y*)A(y*, y), and substituting such 
estimates into (18. 16[) recursively we see that 

f t (x*, y*, y) > Ft\x\ y*)A(y*, y), n = 0, 1, 2, . . . 
On the other hand, we know that 

]imF t H (x*,y*) = P*(t;x*,y*), 

n— >oo 

see Section HJ |Fel40] . |An d91j. Hence, any nonnegative solution of (18. 8p . ( 18. 9 p is 
bounded by P*(t; x*, y*)A(y*, y) from below, and the proof of Proposition 18.31 is 
complete. □ 

The following statement is the analog of Proposition 2.2 in | BF08+| . 

Corollary 8.4. Letfi*(x*) be a probability measure on E* . Fort>0, let (x*(t),x(t)) 
be an E^ -valued random variable with 

Prob{(z*(t),x(t)) = (x*,x)} = £ V*(y*)A(y\y)P^(t;( y *,y),(x\x)). 

(y*,y)eEW 

Then for any time moments < to < t\ < ■ ■ ■ < tj~ < tfc+i < • • • < tk+i, the joint 
distribution of 

(x*(t ),x*(ti), . . . ,x*(t k ),x(t k ),x(tk+i), • • • ,x(t k +i)) 
coincides with the stochastic evolution of fi* under transition matrices 

(P*{to), P*(ti — to), ■ ■ ■ , P*{tk — tk-i),A, P(tk+i — tk), ■ ■ ■ , P(tk+i — tk+1-1)) 
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Proof. In the joint distribution 

v*(y*)Hy*, y)P*{to; (y*, y), K> x o))P*(ti - t ; (x* , x ), (x*, • • • 

■ ■■P*(t k+ i-t k+ i- 1 ; (xlu^Xk+i-i), (x* k+l ,x k+l )) (8.17) 
one uses (18.71) to sum over y, xo, ■ ■ ■ , %k-i and (18. 6p to sum over x* k+1) . . . , x* k+l . □ 

8.3. Construction of multivariate Markov chains. Let Ei, . . . , En be count- 
able sets, Qx, . . . , Qn be matrices of transition rates on these sets, and A^, . . . , A^_ x 
be stochastic links: 

At_i:£ fe X [0,1], A*_ 1 (x^) = l, x6£ t) k = 2,...,N. 

It is also convenient to introduce a formal symbol Aq with Aj( • , • ) = 1. It can be 
viewed as a stochastic link between E\ and a singleton i?o- 

We assume that for each of the matrices Qj, A^_ 1; each row contains only finitely 
many nonzero entries, and that the following commutativity relations are satisfied: 

A *-i(z> u)Q k -i(u, y)=J2 Gk(*> v ) A k-i( v > V)i k = 2,...,N, 

or AjLjQfc-i = QfcA^j in matrix notation. If A^_ 1 (x,?/) = 0, the terms with u = y 
and v — x give no contribution to the sums and thus can be excluded. In that case 
we define (x G E k , y G E k -i) 

A k -i(x,y) := A k k _ 1 (x,u)Q k ^ l {u,y) = Qk{x,v)A k k _ l {v,y), (8.18) 
and also 

( Qk(x,v)A k k _ 1 (v,y) k , 

Q k (x,v,y) = < Al^ix^y) 

[O, if A*_ 1 (a; > y) = 0. 

In case A^_ 1 (x,y) 7^ 0, Q k (x,v,y) is a probability distribution in v G -E* that 
depends on x and y. 

In the application of this formalism that we consider in the next section, there is 
always exactly one v that contributes nontrivially to the right-hand side of (18.181) . 
which means that the distribution Q k (x,v,y) is supported by one point. 

We define the state space E^ for the multivariate Markov chain by (18. lfl and 
then define the off-diagonal entries of the matrix of transition rates on E^ as 
(we use the notation X N = (xi, . . . , xn), Y n = (yi, ■ ■ ■ , Vn)) 

A k-i(yk,x k -i) 



Qk{x k ,y k ) 



I K \^ Kl UK) , k 1 

D( n )(Y V \ — ) I *- k -i\X k ,X k --L / 



Qkixk.ykj^i — 77—: — -Q k +i{x k+ i,y k+ i,y k ) ■ ■ -Qi{xi,yi,yi-i), 
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where for the first line we must have Xj = yj for all j ^ k and some k = 1, . . . , N, 
while for the second line we must have Xj = yj iff j < k or j > I for some 1 < k < 
I < N, and A(xj, yj-i) = for k + 1 < j < I. If neither of the two sets of conditions 
is satisfied, we set Q( n \Xn,Yn) to 0. 

The diagonal entries (Xn, XjA are defined by 

qW(x s ,x n ) = - Q {n) (*n,y n ). 

Y N ^X N 

The definition of can be interpreted as follows: Each of the coordinates x k , 
k = 1, . . . , N, is attempting to jump to y k G E k with certain rates. Only y k s with 
Q(x k ,y k ) 7^ are eligible. Three situations are possible: 

(1) The change of x k to y k does not move X N out of the state space, that is 
Al +1 (x k+1 ,y k )Al_ 1 (y k ,x k ^ 1 ) ^ 0. Such jumps have rates Q k (x kl y k ) i^t^TVT'S - 
Note that for k — 1 the last factor is always 1. 

(2) The change of x k to y k is in conflict with x k -i, that is A-k-iiVk, £fc-i) = 0. 
Such jumps are blocked. 

(3) The change of to y k is in conflict with x k+ i, that is A^ +1 (xfc + i, y^) = 0. 
Then x k+ ± has to be changed too, say to y k +i- We must have A k k +1 {y k+ i, y k ) ^ 0; 
relation (18. 18ft guarantees the existence of at least one such y k+ \. If the double jump 
(x k ,x k+ i) — ¥ (y k ,y k+ i) keeps X N in the state space, it is allowed, and its rate is 

Qk(xk,yk) t^lTT' 1 ] Qk+i(x k+ i } y k+ i,y k ). Otherwise, x k+2 has to be changed as 
well, and so on. 

To say it differently, unless A^_ 1 (|/ fc , x k -\) = 0, the move x k —¥ y k always happens 

with rate Qk{x k ,y k ) jjF^f^-^-^ , an d it ma y cause a sequence of displacements of 

x k+ i, x k+ 2, • • • , where each next Xj uses the distribution Qj(xj, ■ , yj-i) to choose its 
new position. Displacements end once X^ is back in E^ N \ This description implies 
the following formula for the diagonal entries of Q^: 

k=ly k €E k :y k ^ k A h-l&k, Xk-l) 

The definition of is explained by the following statement. 

Proposition 8.5. Consider the matrix A with rows marked by elements of En, 
columns marked by E^' 1 ^, and entries given by 

A(x N , (xx, . . . , x N -i)) = A^_ x (xiv, xn-i) ■ ■ ■ Ai(x 2 , xi). (8.20) 

Then the commutativity relation AQ^ N ~^ = QnA holds. 
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Proof. We have 

AQ^W^-i) = A(x^,X JV _ 1 )g( 7V - 1 )(X w _ 1 ,F Ar _ 1 ). (8.21) 

By (18.191) . the contribution of X^_i = Y^-i to the right-hand side has the form 

N ~ 1 \k ( \ 

fc=l z k &E k :z k ^y k fe ' yfe i; 

For Xjv_! 7^ Vjv-ij the contribution of matrix elements of Q^^pC/v-i, Yjv-i) that 
correspond to jumps . . . ,x t ) -)■ y fc +i, ■ ■ - ,yi), ^- < k < I < N, with all 

other = has the form 

E A^^xjv, yjv-i)A^(yjv-i, 2/JV-2) • • • A]+^(y H2 , yj+i) 

x A} +1 (j/j+i, Xi) A '-i(^i> ^-0 • • • Afe +1 (x fc+ i, x fc ) 
x At-i(zfc, y*-i)A{l^(y fc _i, y fc _ 2 ) • • • A?(jfe, y x ) 

x Q k (x k ,y k ) . k k ' 1< f k ' Vk \ Q k +i{xk+i,yk+i,yk) ■ ■ -Qi(xi,yi,yi-i), (8.23) 

A k-A X k> Vk-l) 

where the summation is over x k , ■ ■ ■ , xi satisfying Xi 7^ y-i for all k < i < I and 

Ati(**> ^ 0, Atx^i, = 0, fc < 1 < I. (8.24) 

Denote this expression by A(k,l). 

Observe that in (|8.23|) . the factors A^_ 1 (xfc, y k -i) cancel out. Let us denote by 
B(k,l) the sum of same expressions (18.231) with canceled A^L^a^, yk-i), and with 
conditions (18.241) replaced by 

^t-iixk, Vk-i) = 0, Aj_ 1 (ar 4 , 2/<_i) = 0, k<i<l. 

Thus, the sum A(k, I) + B(k, I) has no restrictions on x k other that 7^ 
Using the definitions of A^ +1 and Qfc+i we see that 



2J A^ +1 (x fc+ i, x k )Q k (x k , yk)Qk+i(xk+i,yk+i, Vk) 

x k --x k ^y k 

= Qk+i( x k+u yk+i)A- h k +1 (yk+i, Vk)- (8.25) 
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Hence, A(k,l) + B(k,l) = B(k + 1,1). Noting that B(l,l) = 0, we obtain, for any 
Z = l,...,jV-l, 

A(l, I) + A(2, /) + ■•• + A(l, I) = A(l, I) + B(l, I) 
A(x N , Y/v-i) 



A\ +1 {y l+1 , yi ) <^ 



4 5^ A ! +1 (^+i ) ^) ( 5i(^>2/0 - H x N,YN-i)Qi(yi,yi) 



A(x N ,Y N - 1 ) I - Qi(yi,yi) 



^ El+1 Aj+W,v«) 

A l +\z l+l ,yi)_ 

yi) 



+ Q {+ i(yi +1 , y J+1 ) - yi)j , (8.26) 

where we used the commutativity relation Aj +1 Qz = Qi + iA\ +1 along the way. Hence, 
using (I8.2ip we obtain 

aq^W^-i) = A(A ; ,z) + A(^,y JV _ 1 )g( iV - 1 )(iv_ 1 ,yv„ 1 ) 



KJfe<KJV-l 



^ Qn(xn, z n )A(z n , Y N -t) + Qn(xn, x n )A(x n , Y/v-i) = QatA^, Y/ 



(8.27) 
□ 

For any N > 2 and a;jv_i G -Eat_i let us define a matrix Q XN _i_ of transition rates 
on the fiber 

E x N -i = i x N <E -Etv I A 1 §,_ x (x N ,x N - 1 ) ^ 0} 

via 



QxN-ii^N, yjv) — Qn(%n, 2/iv) tiv — 7i — \> yN 7^ X N, 



A A r_ 1 (?/Ar, Xjv- 

(8.28) 

Qx N -l (%Ni x n) — ~ / Qx N -! 

The next statement is analogous to Proposition 2.5 in [B F08+ . 
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Proposition 8.6. Assume that the matrices of transition rates Qi, ■ ■ ■ ,Qn and 
Qxn ■ ■ ■ Qx N _! for any Xj G Ej, j = 1, . . . , N — 1 are regular. Then Q^ 2 \ . . . , 
are also regular. Denote by {P/(£)}i<7<iV and p( N \t) the transition matrices for 
{Qj(t)h<j<N andQW(t). 

Let fix be a probability measure on Ejf, and for t > 0, let (xi(t), . . . ,xjf(t)) be a 
E^ -valued random variable with 

Prob{(a;i(t), . . .,x N (t)) = (xi, . . .,x N )} 

= Yl m(VN)HyN,Y N ^)P^(t;Y N ,X N ). (8.29) 

Y N £EW 

Then for any sequence of time moments 

o<^<t^<...<tl=tr i <t^ 1 <---<tl-_\ = tt 2 <... 

...<tl 2 =tl<t\<---<t l ki (8.30) 

the joint distribution of {x m (t™)} ordered as the time moments coincides with the 
stochastic evolution of under transition matrices 

P N (t ), Pjv(£i -t ),..., PN(tk N ~ h N -l)j^-N-lJ 

p +N-1\ p (+N-1 ,N-1 \ aJV-1 

Pv-i(ti -t ),---,-t J N-i{t kN _ 1 -t kN _ i _ 1 ),A N _ 2 ,... 

...,P 1 (t\-t l ),...,P 1 (tl l -tl l _ 1 ). (8.31) 

Proof. It is a straightforward computation to see that the construction of the bivari- 
ate Markov chain from the previous section applied to Q — Q^ N ~ X \ Q* = Qn, and 
A given by (I8.20p (the needed commutativity is proved in Proposition 18.51) . yields 
exactly Q^ N \ We apply Corollary 18.41 and induction on N concludes the proof. 

□ 

Corollary 8.7. In the assumptions of Proposition \87b\ (P^ N \t))t>o is central in 
the sense of Definition ^. 1\ and the induced semigroup on E N is exactly (.P/v(t))t>o- 
Furthermore, compatibility relations of Proposition \8.S\ also hold. 

Proof. The first two statements follow from Proposition 18.61 with 

k N = 1, k N _i = k N _ 2 — ■ ■ ■ — k\ — 0. 

The third statement is (JS2D with Q = Q^- x \ Q* = Q N , Q {2) = Q (JV) , and A given 
by (jg^Dj) . □ 

9. Stochastic dynamics on paths. Gelfand-Tsetlin graph 

9.1. Central measures on paths and the boundary. Let us return to our con- 
crete setup, cf. Sectional We have E N = GT N , the space of signatures of length 
N, and E^ of (18.11) is the set of Gelfand-Tsetlin schemes of length N; we denote 
it by GT {N} . 
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Due to (HO), the notion of centrality for ^ G M P (GT {N) ) means the following, 

cf. (E2D: For any A = (A (1) -< A (2) -< < A (7V) ) G GT (A °, ^ N \\) depends only 

on A^. For branching graphs, the notion of central measures was introduced in 
[VK81j . see also |Ker03j . 

In Subsection 18 . II we explained that central measures on the space 

E (oo) = . GT (oo) 

of infinite Gelfand-Tsetlin schemes are in bijection, thanks to Theorem I3.1[ with 
Ai p (Q). Let us make this bijection more explicit. 

Given a signature A G GT^v, denote by A + and A~ its positive and negative parts. 
These are two partitions (or Young diagrams) with £(X + ) + ^(A~) < N, where £( ■ ) 
is the number of nonzero rows of a Young diagram. In other words, 

A = (A+ . . . , A+ 0, . . . , 0, -Ar, . . . , -Ar), k = £(X + ) } I = £(A~). 

Given a Young diagram u, denote by d(v) the number of diagonal boxes in v. 
Introduce Frobenius coordinates of v via 

Pi(v) = Vi - i, q i (v) = v' i -i, i = l...,d(u), 

where v ' stands for the transposed diagram. We also set 

Viiy) = Qi(v) =0, i> d{v). 

An element A = (A^ -< A (2) -<...) e GT (oo) , which can be viewed as an infinite 
increasing path in the Gelfand-Tsetlin graph GT, is called regular if there exist 
limits 



a. i = lim — — , & = lim — — , i = l,2,..., 5 = lim 



iV-s-oo N N^oo N N^oo N 

The corresponding point uj = (a ± , f3 ± , 5^) G Q is called the end of this path. 

Theorem 9.1 ( |Ols03j ). Any central measure on GT 1 - 00 -* is supported by the Borel set 
of regular paths. Pushforward of such measures under the map that takes a regular 
path to its end, establishes an isomorphism between the space of central measures on 
GT {oo) andM p {Sl). 

We refer the reader to Section 10 of [Ols03] for details. 

9.2. Matrices of transition rates on GT (A °. With E N = GT N , = GT {N) 
and Qn = V^ N \ let us write out the specialization of the matrix Q from Subsec- 
tion |H2J We will use the notation for the resulting matrix of transition rates 
on GT^. As for the parameters, we will use (16 .ip and (16.21) as before. 

To any A G GT^ we associate an array {l\ | 1 < i < j, 1 < j < N} using 
l\ = + j — i. In these coordinates, the interlacing conditions A^ -< A^ +1 ^ take 
the form 

H +1 >4>iftl 
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for all meaningful values of i and jE 

Similarly, assign GT^ 3 u < — > {n{ = + j — i}i<i<j,i<j<N- Gathering all the 
definitions together, we obtain that the off-diagonal entries of have the form 



DW(A, 



V) 



{lf-z-k + l){lf-z'-k + l), 
(lf + w)(lf + w'), 



where for the first line we must have i, k and l,l<i<k<l<N, such that 
tk-^lHk-j, »1, , l{ +k j • 1 for all k<j<l, 

and all other coordinates of A and v_ are equal, while for the second line we must 
have i, k, I with 1 <i < k < I < N such that 

l{ = if , n \ = l\-l for all k<j<l, 

and all other coordinates of A and v_ are equal. 

The Markov chain generated by can be described as follows: 

(1) Each coordinate If tries to jump to the right by 1 with rate (If — z — k + 1) (If — 
z' — k + 1) and to the left by 1 with rate {if + w)(lf + w'), independently of other 
coordinates. 

(2) If the /f-clock of the right jump rings but If = ifll, the jump is blocked. If its 
left clock rings but If = if" 1 + 1, the jump is also blocked. (If any of the two jumps 
were allowed then the resulting set of coordinates would not have corresponded to 
an element of GT^ as the interlacing conditions would have been violated.) 

(3) If the right if -clock rings and there is no blocking, we find the greatest number 
I > k such that l\ = if + k — j for j = k, k + 1, . . . , I, and move all the coordinates 
{l{} l j =k to the right by one. Given the change If h- > If + 1, this is the minimal 
modification of the set of coordinates that preserves interlacing. 

(4) If the left If -clock rings and there is no blocking, we find the greatest number 
/ > k such that l 3 i+ j_ k = If for j = k, k + 1, . . . , I, and move all the coordinates 

{H+j-kVj=k t° the left by one. Again, given the change If !->■ If — 1, this is the 
minimal modification of the set of coordinates that preserves interlacing. 

Certain Markov chain on interlacing arrays with a similar block-push mechanism 
have been studied in |BF08+| , see also |BK10j . In those examples the jump rates 
are constant though. 



One could make the interlacing condition more symmetric (both inequalities being strict) by 
considering the coordinates l\ = + (j + l)/2 — i instead. This would imply however that 
lj G Z + 1/2 for odd j while i| S Z for even j. 
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9.3. Regularity. In order to claim the benefits of Proposition 18.61 and Corollary 
18.71 we need to verify the regularity of the fiber matrices of transition rates (18.281) . 
In our concrete realization, they take the following form. 

For any N > 2 and any k G E^-\ — GUV-i, the fiber E K =: GT K C GT/v takes 
the form 

GT K = {A G GT N | k -< A}. 

Using the coordinates {li = N + Aj — i}f =1 for A G GT N and {rij = AT + ^ — 

for z/ G GTjy, the off-diagonal part of the matrix of transition rates T> K := Q K on 

the fiber GT„ has the form 



Proposition 9.2. For any N > 2 and any k G GTjy-i, i/ie matrix of transition 
rates V K on GT K is regular. 

Proof. The interlacing condition in the definition of GT K implies that T> K is the 
matrix of transition rates for independent birth and death processes conditioned 
to stay within A^ non-overlapping intervals inside Z; one interval per process. The 
results of Section 3.2 show that any such birth and death process is regular as such 
a process either lives on a finite set or it is a one-sided birth and death process of 
the type considered in the proof of Theorem 15.11 □ 

Corollary 9.3. For any N > 1, the matrix of transition rates on GT^ 

is regular, and the corresponding semigroup (P'- 7V ^(t))t>o is central. The induced 
Markov semigroup on G1V coincides with that of Section^ 

Proof. Follows from Proposition 18.61 and Corollary 18.71 □ 

9.4. Exclusion process. Observe that the projection of the Markov chain gener- 
ated by to the coordinate l\ is a bilateral birth and death process. Further- 
more, the jumps of l\ are only influenced by l\, the jumps of If are only influenced 
by l\ and l\ , and so on. On the other side, the jumps of l\ are only influenced by 



{llJl.. .l k k Z\} for any k> 2. 

Hence, the projection of the Markov chain defined by to the coordinates 

(In < l N-l < • ■ ■ < l \ < l l < ■ ■ ■ < l i) is also a Markov chainE The fibers of 




where for the first line we must have i, 1 < % < N, such that 

Hi = li + 1, rij = lj for j 7^ i, 



and for the second line we must have 

Hi = U — 1, nj = lj for j ^ i. 



2 Once again, all the inequalities would be strict if we considered coordinates l\ — \\ + (j+l)/2— i. 
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this projection are finite, hence, according to Proposition 14.44 our Markov chain on 
GT^- 1 collapses to the smaller one, whose matrix of transition rates is also regular. 

Let us project even further to (l\ < l\ < ■ ■ ■ < Killing extra coordinates one 
by one and using the results of Section loTTl to verify the regularity for the fiber chains, 
we see that the collapsibility of Proposition 14.41 holds. Let us give an independent 
description of the resulting Markov chain on {l{}j>i- 

Set 

2)iv = {yi < V2 < ■ • ■ < Vn I Vj e Z, l < 3 < N}, 
2)oo = {yi < y 2 < ■ ■ ■ I yj e Z, j > l}. 

Define the matrix D^l of transition rates on 2)^ by 

D (w)/ y / y,n = Uy k -z-k + l)(y k - z' - k + 1), 
\(yk + w)(y k + w'), 

where for the first line we must have k and /, 1 < k < I < N, such that 

y'j=y'k + k- J, y'- =y'j + l for all k<j<l, 

and all other coordinates of Y' and Y" are equal, while for the second line we must 
have 

y'L = yi - !> y'L = yL fo r m + k - 

In other words, each coordinate tries to jump to the right by 1 with rate 
(y k — z — j + l)(yk — z' — j + 1), and it tries to jump to the left by 1 with rate 
{yk + w){yk + w'), independently of other coordinates. If the left y k clock rings but 
yu = yk-i + 1 then the jump is blocked. If the right 2/fc-clock rings we find the 
greatest number / > k such that yj = y k + k — j for j = k, k + 1, . . . , I, and move all 
the coordinates {y k , ■ ■ ■ ,yi} to the right by one. One could think of y k "pushing" 
y k +i, ■ ■ ■ ,yi- Alternatively, if one forgets about the labeling one could think of yk 
jumping to the first available site on its right. 

Clearly, these Markov chains are compatible with projections 2)jv+i — > 2) at that 
remove the last coordinate. Thus, we obtain a Markov semigroup on h_m2)Ar = 2)oo. 

This semigroup is a sort of an exclusion process — it is a one-dimensional interact- 
ing particle system with each site occupied by no more than one particle (exclusion 
constraint). A similar system, but with constant jump rates, was considered in 
[BF08] and called PushASEP. A system with one-sided jumps and blocking mecha- 
nism as above is usually referred to as Totally Asymmetric Simple Exclusion Process 
(TASEP), while a system with one-sided jumps and pushing mechanism as above 
is sometimes called long range TASEP. See |Spi70| , |Lig99] for more information on 
exclusion processes. 

Proposition 9.4. The exclusion process defined above has a unique invariant prob- 
ability measure. With probability 1 with respect to this measure there exists a limit 
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r = linijv^oo un/N, which is a random variable with values > 1. Under certain 
additional restrictions on parameters (z, z' ,w,w') , see below, the function 

a(s) = s(s — 1)— Prob{r < s} — a\s + |(d3a 4 + a\) 
is the unique solution of the {2nd order nonlinear) differential equation 

— a'(s(s — l)cr") = (2((s — |)cr' — a)a' — a^a^a^ 2 

- (a' + a?) (a' + + a 3 ) V + a\) (9.1) 

with boundary condition 

a(s) = -a\s + ^(a 3 a 4 + a{) + = s" 2ai + o(s iai ), s ->■ +oo, 

where the constants a±, 02,03,04 are gwen fry 

Z + z' + 10 + If' Z — + U> — w' 2 — 2' — ttf + «/ 

ai = a 2 = , a 3 = , a 4 = . 

Remarks 1. The quantity lini7v->oo Vn/N can be viewed as the asymptotic density 
of the system of particles {yf) at infinity. Proposition [931 claims that for the invariant 
measure, this quantity is well-defined and random. 

2. The restrictions on parameters come from Theorem 7.1 of |BD02j . They can 
be relaxed, see Remark 7.2 in [BD02] and the end of §3 in | Lis09+ . 



3. The differential equation above is the so-called cr-form of the Painleve VI 
equation first appeared in [J M81] . 



Proof of Proposition 9.4 ■ The invariant measure is simply the projection to 

yi = A«, y 2 = \f ) + l, y 3 = \f + 2, ... 

of the central measure on GT*- 00 -* corresponding to the spectral zw-me&sme. The 
uniqueness follows from the uniqueness of invariant measure on countable sets 2) n, 
cf. Theorem 1.6 of |And91] (a similar argument was used in the proof of Theorem 
I7.2p . The existence of liniN^oo t/Ar/iV follows from Theorem 19. II Finally, the charac- 
terization of the distribution of this limit in terms of the Painleve VI equation was 
proved in Theorem 7.1 of |BD02j . see Lis 09+ | for another proof. □ 

10. Appendix 

10.1. Truncated Gelfand-Tsetlin graph. Fix two numbers k, I — 0, 1, 2, . . . not 

equal to simultaneously. Denote by GT N (k, I) the subset of GT N formed by the 
signatures A subjected to the restrictions 

k > Ai > •• • > Ajv > -I. 

Obviously, this subset is finite and nonempty, and if A G G1V(&, I) and v -< A then 
v G GT/v-i(&, I). Thus, the union of the sets GT N (k,l) for N = 1,2,... forms 
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a subgraph of the Gelfand-Tsetlin graph. Let us denote this truncated graph as 
GT(fc, I). The definition of the links A^ +1 (A, v) in the truncated graph remains the 
same, the only difference is that we assume A and v to be vertices of GT(k, I). 
The boundary of GT(k, I) is the subset Q(k, I) C Q determined by the restrictions 

a ± = 0> 7 ± = 0, 0+ = O (i>k), /3r=0 (j > /), 

so that only k + I parameters are nontrivial: 

0t > ■ ■ ■ > 0t > 0, 0i > ■ ■ ■ > 0i > 0, 

where, as before, 0f + 0i — 1- The definition of the links A^(w,A) remains the 
same, only u is assumed to belong to Q(k,l). Viewing Q(k,l) as a subset of M. k+l 
one sees that it is a closed simplex of full dimension. 

The boundary is Feller. Since Q(k, I) is compact, this simply means that the links 
A°°((X>, A) are continuous in w 6 Q(k, I). 

Fix parameters 

z = k, z' = k + a, w — I, w' — l + b, where a, b>— 1. (10.1) 

Setting 

u = z, u = z', v = w, v' = w' 

in (15. ip we obtain a truncated birth and death process on {— /, . . . , k} C Z; the rates 
of jumps k — > k + 1 and — / — )■ — / — 1 being equal to 0. Note that the nonnegativity 
of the jump rates on {— /, ...,&} C Z is ensured by the inequalities a, b > —1. 

More generally, for any N > 1 the same expressions as before correctly determine 
a matrix of transition rates on GjTjv(/c, I). Due to finiteness of the state space, the 
existence of the corresponding Markov semigroup (-P/v(t))t>o becomes obvious. 

The semigroups (Pjv(0)t>o with varying N = 1,2,... are consistent with the links 
and thus determine a Feller Markov semigroup (P(t)) t >o on the boundary Q(k,l). 

The expression for M Z)J5 / )11 , j . u ,/|jy(A) given in §5.1 vanishes unless A belongs to the 
subset GTjv(fc, J), and it is strictly positive on this subset. Thus, the same defi- 
nition gives us a probability measure on GTjv(fc, Z). The invariance property and 
the compatibility with the links A^ +1 remain valid. The limit measure lives on the 
boundary Q(k,l) and it is invariant with respect to the semigroup (P(t)) t >o. 

The whole picture sketched above is consistent with the automorphism of GT 
described in Remark l3.21 More precisely, the shift of all coordinates of signatures by 
1 amounts to the transformation k — > k + 1, I — > I — lof the main parameters (the 
other two parameters a and b do not change). The use of this shift automorphism 
allows one to reduce the case of general parameters (k, I) to the special case (n, 0) 
with n = k + I, which simplifies some formulas and computations. 

Hence, let us now assume that / = 0. Then the coordinates 0~ disappear and we 
are left with n coordinates yi := (3f subjected to 

1 > Vl > ■ ■ ■ > Vn > 0. 
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The signatures A G GTn(u, 0) may be identified with Young diagrams contained 
in the rectangular shape Nxn; that is, A has at most N rows and n columns. Under 
this identification, one has a simple expression for the link: 

N , x 

AS(w, A) = Dim^ A TT(1- y t ) N sy . . . , (10.2) 

ti V 1 — 2/1 l-J/n/ 

where s\> is the Schur polynomial indexed by the transposed diagram A'. 
Further, the invariant measure on the boundary takes the form 

n 

const- J] (yi-ytf-llil-yiyytdyi (10.3) 

l<i<j<n i=l 

with an appropriate constant prefactor that turns the measure into a probability 
distribution. The random n-tuple (yi, . . . , y n ) C [0, 1] with this distribution is known 
under the name of the n-particle Jacobi orthogonal polynomial ensemble. 

The fact that the integral of (110.21) against the distribution (110. 3p reproduces the 
measure M nin+aj0j b|jv on GT N (n, 0) can be verified directly; this is a version of the 
Selberg integral. 

For more detail, see [Ker03j and |BO05aj . 

Since the boundary has finite dimension, there is a possibility to describe the 
Markov process defined by (P(t)) t >o more directly; this is done in the theorem 
below. 

Consider the ordinary differential operator associated with the Jacobi orthogonal 
polynomials with weight (1 — y) a y b , 

= y (l-y)^ + [b + l-(a + b + 2)y)^. 

More generally, abbreviate 

V n = V n (yi,...,y n ) = ] [ (yi-yj) 

l<i<j<n 

and consider the partial differential operator in variables yi, . . . , y n given by 

D t b} : =v n °(£; - Vi) Wt + [b + 1 " (a + b + 2)Vi] w) ) ° v " + <: " ' 



b + 1 - (a + b + 2)yt + 



2yi(l - yi 



Vi - Hj 



d_ 

dyi 



where (• ■ ■ ) stands for the constant annihilating the constant term arising from the 
conjugation by the Vandermonde determinant: 

n-l 

(■ • • ) = m{yn + a + b + 1). 



m=0 
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Although the coefficients in front of the first order derivatives have singularities on 
the hyperplanes y^ = yj, the operator is well defined on smooth symmetric functions 
in variables yi, . . . ,y n , and it preserves this space. It also preserves the space of 
symmetric polynomials. 

Theorem 10.1. Let z = n,z' = n + a, w = 0, w' = b, where a,b > — 1. Let 
yi = 0i , . . . , y n = f3£ be the coordinates on the boundary Q(n,0) of the truncated 
graph GT(n, 0), and recall that 1 > y\ > ■ ■ ■ > y n > 0. 

The Markov process on the simplex Q(n, 0) determined by the Markov semigroup 
(P(t)) t >Q is a diffusion whose infinitesimal generator is the differential operator 

D^' b ^ with an appropriate domain containing the space of all symmetric polynomials 
in variables yi, . . . ,y n . 

As shown in |Gor09] . the same diffusion process also arises in a scaling limit 
transition from some discrete time Markov chains on the sets GIV(n, 0), as iV — > oo. 

The above description of the infinitesimal generator remains true in the general 
case (110. ip of the degenerate series parameters. The only change concerns the cor- 
respondence between the /^-coordinates and the ^/-coordinates; now it takes the 
form 

( Vl , . . . , y n ) = (1 - 0f, . . . , 1 - /3f, 0+, . . . , P+), n:=k + I. 

10.2. The formal generator. A natural question is how to extend the explicit de- 
scription of the infinitesimal generator obtained in Theorem IIP. II to the case of gen- 
eral (admissible) values of parameters (z, z', w, w'). There are some indications that 
a direct generalization is impossible, in the sense that the generator cannot be ex- 
pressed as a second order differential operator in the natural coordinates (cr^ 5^) 
on Instead of this, we present below an explicit expression for the generator in 
a different system of coordinates. 

Introduce the functions (p n = ip n {u) on fl, n 6 Z, as the coefficients in the Laurent 
expansion of ^^(u), see (13.31) and (13. 4p . 

The functions ip n are continuous and nonnegative on Q, and they satisfy the relation 

y^y n (q;) = 1, u E Q, 

which is an immediate consequence of the fact that $ w (m) takes value 1 at u — 1 for 
all uefi. 

Let us extend this definition by setting, for any iV = 1,2,... and any signature 
v E GT N , 

(p v (u) = det[ Vui _ i+j {u)]. J=1 , u e Q. 



3 A possible explanation is that these coordinate functions are not in the domain of the generator. 
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The functions <p u appear in the expansion 

N 



i=X 



Their fundamental role is explained by ( 13. 2p . 

By definition, the functions ip v are contained in the algebra of functions generated 
by the functions tp n , iifZ. Conversely, any monomial in ip n 's of degree N can be 
expanded into a series on the functions ip v , v G GIV- Namely, for arbitrary integers 
rii, ... , tin, one has 

<Pni---<Pn ir = ^2 K i v \ n li - ■ - i n N)¥v, 

where the numbers K(. . . ) are defined as the "rational" analogs of the Kostka num- 
bers, that is, these are the coefficients in the Laurent expansion of the rational Schur 
functions, 

s„(ui, . . . ,u N ) = 22 K(iy\n 1 ,...,n N )u r l 1 ...u N N . 

ni,...,njv€Z 

For instance, for N = 2 and v — (vx, v 2 ) we have 



K(u x , v 2 I nx,n 2 ) 
which implies 



1, if rii + n 2 = v\ + v 2 and Ini — n 2 \ < v\ — u 2 , 



0, otherwise, 



L Pn 1 i Pn 2 — ^ ] l -P(n 1 +p,n 2 -p)- 
p=0 



Definition 10.2. Fix an arbitrary quadruple (z, z' , w, w') of complex parameters 
and introduce the following formal differential operator in countably many variables 
{tp n : n G Z} 



ni>n2 



'd(p ni d(p ri 
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where, for any indices n\ > n 2 , 

oo 

= ^( n i ~n 2 + 2p + l)(ip ni+p+1 tp n2 _ p + V9 ni+p v5„. 2 -p-i) 

p=0 

oo 

-(ni - n 2 )(Pni ( Pn2 - 2 X^™ 1 ~ H 2 + ' 2 P)Vn 1 +p i Pn2-p 

p=l 

and, for any n G Z, 

-B n = (n + w + l)(n + w' + l)y n +i + (n — z - l)(n — z' - 1 1 
— [in — z)(n — z') + {n + w){n + w'))ip n . 

Note that only coefficients B n depend on the parameters (z,z',w,w'). 

Assume now that the quadruple (z, z', w, w') satisfies the condition (16.11) . As 
above, let P/v(t) (N = 1, 2, ... ) and P{t) be the corresponding Markov semigroups 
and let An and A stand for their infinitesimal generators. These are densely defined 
operators in the Banach spaces C (GT N ) and C (Q), respectively. We know that 
all finitely supported functions on G1V belong to the domain of An- This implies, 
cf. (13.21) . that all the functions ip v lie in the domain of A. 

Theorem 10.3. Let A be the infinitesimal generator of the Markov semigroup 
(P(t)) t >Q with parameters (z,z',w,w') satisfying condition (16. ip . and ID) be the for- 
mal differential operator introduced in Definition \10.2 . Regard D as an operator 
from the space <C[(p n : n G Z] of polynomials in countably many variables (p n , n e Z, 
to the larger space C[[(f n : n G Z]] of formal series in the same variables. 
Then for any N = 1, 2, . . . and any v G GTn one has 

H(p„ = A(p u . 

Moreover, D is the only formal second order differential operator in variables tp n , 
jiGZ, with such a property. 

Note some properties of D: 

1. Formal application of D to the infinite series Ylnez gives 0. This agrees with 
the fact that the sum of this series on Q equals 1 and the fact that Al = 0. 

2. For any fixed integer m, D is invariant under the change of variables (p n — > (p n + m 
[n G Z) combined with the shift of parameters 

z — > z + m, z' — > z' + m, w -» w — m, w' — >■ w — m, 

cf. Remark 3.7 in [BO05aj . 

3. Set z = k and w = I, where (fc, I) is a couple of nonnegative integers not equal 
to (0,0). Then D respects the relations 

• • • = tP-l-2 = <f-l-l = = ip k+1 = ip k+2 = ■ ■ ■ 
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and thus can be reduced to an operator in the polynomial algebra C [ip~i, ■ ■ ■ , fk\- 
More precisely, this means the following assertion: When © is applied to a monomial 
containing at least one variable ip n with index n outside [— I, k] then all monomi- 
als entering the resulting series with nonzero coefficients have the same property. 
Indeed, this follows from the structure of the coefficients A nin2 and B n . 

Moreover, the resulting operator in the algebra C[(p-i, . . . , ip^] can be further re- 
duced modulo the relation 



and then it coincides with the differential operator from Theorem 110. 1[ where 

a — z' — k, b — w' — I. Here we use the fact that quotient algebra 



C[y?_;, . . . , (p h ]/ H h <fk = 1) 

can be identified with the algebra of polynomial functions on the simplex Q(k, I). 
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